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Abstract 

The equations for the electromagnetic field in an anisotropic media are written in 
a form containing only the transverse field components relative to a half plane bound- 
ary. The operator corresponding to this formulation is the electromagnetic system's 
matrix. A constructive proof of the existence of directional wave-field decomposition 
with respect to the normal of the boundary is presented. 

In the process of defining the wave-field decomposition (wave-splitting), the re- 
solvent set of the time-Laplace representation of the system's matrix is analyzed. 
This set is shown to contain a strip around the imaginary axis. We construct a split- 
ting matrix as a Dunford-Taylor type integral over the resolvent of the unbounded 
operator defined by the electromagnetic system's matrix. The splitting matrix com- 
mutes with the system's matrix and the decomposition is obtained via a generalized 
eigenvalue-eigenvector procedure. The decomposition is expressed in terms of com- 
ponents of the splitting matrix. The constructive solution to the question on the 
existence of a decomposition also generates an impedance mapping solution to an 
algebraic Riccati operator equation. This solution is the electromagnetic generaliza- 
tion in an anisotropic media of a Dirichlet-to-Neumann map. 

Keywords directional wave-field decomposition, wave-splitting, anisotropy, electro- 
magnetic system's matrix, generalized eigenvalue problem, algebraic Riccati operator 
equation, generalized vertical wave number. 



1 Introduction 



Wave field decomposition is a tool for analyzing and computing waves in a configuration 
characterized by a certain directionality. The wave-field decomposition, or wave splitting, 
has been used to separate the wave field constituents which are of importance for the 
analysis on a boundary, both for direct and inverse scattering problems [25, 6, 14, 11,9, 
35, 39] and for the analysis of boundary conditions (see for example [5, 27, 2, 23]). 

A remaining challenge in seismic prospecting methods is to incorporate anisotropy into 
the analysis. The enormous data sets used for studying such inverse problems are on the 
border or beyond today's computers [3] . A common method to access such problems is to 
use wave field approximations. Such approximations have been developed for and applied 
to a wide range of hyperbolic equations describing wave propopagation in isotropic media. 
One class of such approximations is based on a decomposition of the wave-field into up- 
/down-going components. Such a decomposition is usually denoted a wave splitting or 
a wave-field decomposition. There are essentially two types of limitations to the present 
theory of wave splitting: The traditional operator based approach of wave-splitting has 
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been limited to heterogeneous isotropic materials sec e.g., [38, 14, 9, 28] or up-/down 
symmetric media [16]. This method is based essentially on constructing a certain square 
root operator and it fails, once the media becomes inherently anisotropic. Whereas wave- 
splitting by spectral decomposition of a certain matrix is restricted to homogeneous or 
depth-independent material. Here both anisotropic and bi-anisotropic materials have been 
considered [12, 30]. 

The present paper removes both of these limitations. We present a derivation of a 
three dimensional wave splitting for electromagnetic fields in the presence of inherently 
anisotropic loss- less heterogeneous constitutive relations. We show the existence of a 
decomposition by a constructive argument. The decomposition is given for media with 
anisotropic permittivity and permeability that is described by self-adjoint, heterogeneous, 
positive definite matrices. These conditions are sufficient but not necessary material 
conditions for the resolvent set of a certain operator to contain the strip around the 
imaginary axis. The requirements of the material parameters, in time-Laplace domain, 
are corresponding to a medium with only instantaneous lossless response. The analysis 
which use pseudo differential calculus is straightforward if one assumes that the material 
coefficients depend smoothly on the spatial variables. In order to simplify the analysis 
this assumption is made. Physically this should not be regarded as a restriction since the 
smooth functions densely approximate the square integrable ones. 

The decomposition is constructed through a generalized eigenvalue-eigenvector proce- 
dure and a certain commutation of two operators. The construction of the commuting 
operator, the splitting matrix, is made by means of a functional analysis approach using 
the resolvent of the electromagnetic system's matrix and is analyzed with pseudodifferen- 
tial calculus with parameters, to prove its existence and to study its behavior. The method 
is a generalization of the stratified- media case first presented in [12] and extends the the- 
ory form the linear acoustic [22] case to the considerably more complex electromagnetic 
case. The challenge in going from the anisotropic-acoustic case to the electromagnetic 
case includes a more complex differential operator with a non-trivial null-space, as well as 
the analysis of a resolvent operator which here is the inverse of a 4x4-matrix of operators. 

There is a wealth of literature on wave splitting and their applications. A few references 
are mentioned below. For time-domain wave-splitting see [11], where both the wave 
equation and the Maxwell's equations are considered with both applications and theory. 
Wave-splitting in connection with Bremmer series for linear acoustics [14, 8, 9] and uniform 
asymptotics and normal modes [7, 15] has been used to analyze the wave-field constituents. 
An extension to include dispersion is presented in [28] and wave-splitting on structural 
elements in [21]. The square-root of a certain operator is a key step in isotropic wave- 
splitting this operator has been carefully studied in [17]. A reciprocity theorem approach 
to decomposition is used by [33]. The results for anisotropic media includes [12, 16, 22]. 

Applications of the existing wave-splitting techniques include several successful ana- 
lyzing tools of the wave- field including Bremmer series, normal modes and uniform asymp- 
totics. Another spin off is the development of fast numerical codes to calculate the wave 
fields. Among their implementations we have 'rational approximations' and 'generalized 
screens' and 'multiple-forescattering-single-backscattering approximation' [36, 18, 31, 39]. 
In the active field of time-reversal mirrors see e.g., [34, 40], the wave-splitting techniques 
have been used see e.g., [23]. It is our hope that the extension of the wave-splitting 
techniques to the inherently anisotropic case will provide a base for generalization of the 
above mentioned applications to analysis and fast numerical codes to general anisotropic 
media. 

The present paper is organized in a set of three propositions that step by step in- 
troduce and prove the necessary properties and tools to obtain the decomposition. The 
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analysis is preformed in the time-Laplacc domain and the procedures impose Hmitations 
on the Laplace parameter. In §2 the problem is formulated after a rewriting of Maxwell's 
equations to a suitable form. In §3 the properties, mostly the spectral properties, of the 
electromagnetic system's matrix are stated and proved using functional analysis. The 
propositions impose only the natural condition that the Laplace parameter has to belong 
to the right-hand half plane of the complex space. In §4 the splitting matrix is constructed 
and several of its properties are shown. The analysis utilize that the material parame- 
ters are self-adjoint, positive, and furthermore, require a mild constraint on the Laplace 
parameter, in order to obtain a certain ellipticity condition that is needed in the subse- 
quent analysis. The most important property shown in this section is that the generalized 
eigenvectors of the splitting matrix can be obtained explicitly in terms of the elements 
of the splitting matrix. In §5 the decomposition is derived in terms of the generalized 
eigenvectors of the splitting matrix. The last section concludes with a discussion and 
some observations. 

Some lengthy intermediate derivations of the electromagnetic system's matrix in §2 
arc detailed in Appendix A. In Appendix B the special case of an isotropic homogeneous 
medium is treated using the approach developed in the present paper and the results are 
compared with traditional methods. In Appendix C the determinant of the symbol of 
the electromagnetic systems matrix is given. Furthermore, explicit integrations of the 
resolvent in symbol representation are presented in terms of residue calculus and the 
integrals are stated in terms of the roots of the determinant of the principal symbol of 
the electromagnetic system's matrix. 



2 Directional wave-field decomposition 
2.1 The two-way equations for Mctxwell equations 

We consider electromagnetic wave motion in heterogeneous anisotropic media with in- 
stantaneous response. Let a; e be a point in space and f e M is time. The media is 
assumed to be independent of time. The following initial conditions of the fields ensures 
causality, 

E{x, t)=0 , D{x, t) = for f < and all a; e , (2.1) 
H{x, t) = , B{x, t) = for t < and all a; e . 

where B = magnetic flux density [T], E = electric field strength [V/m], D = electric flux 
density [C/m^] and H = magnetic field strength [A/m]. The quantities E,D,H,B are 
all functions of space, a; G and time, t G R, with values in R^, we return to which 
function spaces that they belong at a latter point in the present paper. Above we have 
used standard units and notation see [19]. 

The electromagnetic field satisfies the first-order hyperbolic system of partial differ- 
ential equations in time domain. Maxwell's equations see e.g., [13, 20]. We consider 
Maxwell's equations in time-Laplace domain. That is, 

sB + 'V X E = K" , 

(2.2) 

-sD + Vx H = J^ . ^ ^ 

where J*^ = external electric current density [A/m^] and K" = external magnetic current 

density [V/m^]. The external currents are applied, prescribed, sources. The causality of 
the field is taken into account by requiring that all field quantities are bounded functions of 
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the time-Laplacc parameter s, that is in general complex valued and lies in the right-hand 
plane Re {s} > 0. With the specified initial condition (2.1), we have dt s. In this paper 
X = {xi,X2,X3} are right-handed orthogonal Cartesian coordinates. All the subsequent 
analysis is carried out in the domain of {x, s) hence there is no need to distinguish between 
the time dependent field, E(x,t), and the Laplace-parameter dependent field, E{x,s). 

To explicitly introduce the material parameters into the equations we assume the 
following constitutive relations 

B{x,s) = n{x)iinH{x,s) ^^^^ 
D{x, s) = e{x)eoE{x, s) 

where /i = relative anisotropic permeability tensor, e = relative anisotropic permittivity 
tensor, fj,Q = empty space permeability [H/m] and cq = empty space permittivity [F/m]. 
The relative permeability and permittivity are assumed to be self adjoint and positive 
definite 3x3 tensors of second rank, that is, the media under consideration has only 
instantaneous response. Inserting the constitutive relations into the Maxwell equations 
gives 

s^i^ioH + V X E = K\ 

-seeoE + \7 x H = J" . ^ ' ' 

Before proceeding we re-scale and change dimension of the equations analogous to e.g., 
[9, p. 10] to simplify the subsequent analysis: 

Co ^ := iio^o , s := — , H := ^/JI^H , E := ^/e^E , 
Co 

r:=^r, K'-y/^K^ (2.5) 

where cq is the speed of light in vacuum. Upon substituting (2.5) into (2.4), 

siiH + y X E ^ , 

, , (2.6) 
-stE + VxH=r . 

All the following considerations refer to this transformed space and for notational sim- 
plicity we remove the •', but remember the change in dimension, in particular that s 
has dimension m^^, H and E have dimension (J/(m'^))^/^, J° and K"^ has dimension 
(J/(m )Y^^. The transformation above is for dimensional convenience, in particular in 
the calculus of pseudodifl[erential operators see §4. 

The 'evolution' of the wave field in space, along a direction of preference, can be 
expressed in terms of the change of the wave field in the directions perpendicular to it. 
The direction of preference is taken to be along the a;3-axis (or 'vertical' axis) and the 
remaining ('horizontal') coordinates are denoted by x^, x^, a, f3 G {1,2} or x' = {xi,X2} 
when convenient. The procedure requires a separate treatment of the vertical components 
of E and H. From (2.6) we find the vertical field components to be 

s^aaHs = -sna^H^ - (V x E)^ + , 

sesaEs = sespEf) + {V x H)^ - , ^ ' ' 

where Einstein's summation convention for Cartesian tensors are employed for repeated 
indices a, (3 e {1, 2}, e.g., ^i^pHp = Y%=i /"3/3^/3- To project out the third component of 
a vector we have used the subscript 3, to explicitly show the notation consider 

(V X E)^ := diE2 - d2Ei . (2.8) 
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Thus, Eq. (2.7) relates the vertical components of the electric and magnetic field strength 
to the horizontal components. The remaining equations contained in (2.6) are 

sn^jjHp + sHa^H^i + (V X E)o, = Kl, 
-se^fjEis - se^sEs + (V x iJ)„ = . 

By replacing Es, H3 in (2.9) with (2.7) we arrive (the derivation is detailed in Appendix 

A) to 

{Id3 + A)F = N, (2.10) 

where AT is a linear combinations of the sources and their derivatives cf. (A. 12) and where 
the elements of the electromagnetic field matrix, F, are given by 

Fi := El , F2 := -E2 and F3 := H2 , F4 := Hi . (2.11) 

To simplify some of the following calculations we introduce the notation 

E:={Ei,-E2f and H:={H2,Hif . (2.12) 

There are several possible orderings of the transverse components of E, H in F. The par- 
ticular choice of combinations given in (2.11) has two advantages. First, the given choice 
ensures that both the matrix operators A12 and A21 are invertible. Secondly, we have that 

the third component of the Poynting vector, equals {E x _ff)3 = E'^H = {Fi , -F2) (^3 , -^4)"^, 
where ~ denotes the complex conjugate. The 4x4 electromagnetic system's matrix. A, is 
here represented by four 2x2 blocks 



where each block-matrix is given by 



1 f H23di \ . f dlCsi -81632 \ -I 



■^11 - M33 I ^^3^2 ^13^1 i + I -^2631 ^2632 j > 



in which 



t = s I ^22 T^2i\ _ ^-1 f die^^di -91633^92 
' 2^12 vii ) V -'92633^(91 026^^82 



^21 = s ( I - s^M ^^^^^ ^^'^^^ 

-£21 £22 / V 9iiJb^^d2 din^^di 



^ ^ I ^2^*32 92M31 \ -1 _|_ ^-1 f £1381 -£13^2 
^ l9lM32 81II31 J '^^^ I _g235i £23192 



Saf3 — ^af) — ^a3^33^3f) j 
l^al3 = Ma/3 — IJ'a3 1^33 1^3(3 ■ 



(2.14) 



(2.15) 



The permeability and permittivity are symmetric, 3x3 tensors of rank (tensor order) 2 
that are bounded from below and from above. Hence the upper-left 2x2 matrices of e 
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and ^ are bounded below by the constants eo and jlo respective. To show the notation 
we have 

Uaeal3Ui3 > CoUaUa , (2.16) 

for any complex field u = (ui, U2)- Prom the definitions of Saf} and it is clear that they 
arc symmetric matrices (since e and arc symmetric and all elements are real valued. 
Furthermore, each is bounded below by the constants ei and /<i respectively. This follows 
from the identity (summation over repeated index a,(3 £ {1, 2} and j, k G {1, 2, 3}) 

VjCjkVk = UaeapU/S , (2-17) 

where 

Vi = Ua{S^a - eJa^esaJsi) , for i = {1, 2, 3} , (2.18) 

for any complex field u. Since e is positive definite s must also be positive definite and 
analogously for fj, and v. 



2.2 Preliminaries 

We consider the electromagnetic system's matrix and other operators on Sobolev spaces. 
Let H''(M^;C) be the set of functions belonging to the Sobolev space of order r e N 
with domain in and values in C, with a weighted inner product to compensate for 
the dimension of the derivative. To extend this scalar space to vectors we introduce the 
notation 

W :=W{R^;C^) (2.19) 

for a 4 X 1 matrix, with each element in the set of functions belonging to the Sobolev 
space of order r. Let F = {E'^, H^) and G = {E^, H^). Then we define the inner product 
to be 

{F,G\ = j d?x' J2 (d'.'Ef d',,E^ + d',,EI dl,El 

where z/q is a constant of dimension length and it is used to normalize the change of 
dimension from the derivatives. All components of the field depend on x, but inner 
products and norms refer to x' and we treat 2:3 as a parameter. We have adopted the 
multi-index notation of pseudodifferential calculus [32] above and use k & together 
with 

|fc| = fci + fc2 . (2.20) 
The norm corresponding to the inner product is 

\\F\\r=^/{F,F)r . (2.21) 

The set H'' with the inner product (•, ■)r is a Hilbert space. For the case r = we recover 
the Lebesgue space of square integrable functions, {L^, (•, •)o}. 

The analysis of unbounded operators such as the electromagnetic system's matrix 
— requires that one specifies the domain of the operator and its embedding space. Below 
we consider operators on the space {L^,(-,-)o}, that is A : {L^,(-,-)o} — > {L^,(-,-)o}. 
Since the operator is an unbounded operator, we also need to specify its domain, which 
is H^. The domain of the operator on a space is fundamental for the analysis. Here all 
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operators have dense domains and when necessary, with restrictions to dense subsets of 
the their domains for the operation under consideration to be defined. In the case of such 
a restriction we use the notation A\g for the operator restricted to this dense subset of its 
domain and indicate by q what dense subset is understood to be the restricted domain. 

One can also consider As,x as an operator on {H'', {■,-)r}, (if r > this is a restriction 
of the operator defined above) with domain H''"^, and the analysis extends trivially to 
this case. An alternative method was detailed in [22], where one instead of As.x consider 
the operator 'y~^As,x^^, where Y '■= (1 ~ Uo'^aday^^- The results in this paper hold also 
for this class of operators. 



2.3 Formulation of the problem 

To be able to solve the scattering process along the vertical direction separately from 
the scattering process in the horizontal directions, we diagonalize the operator on the 
left-hand side of (2.10). This procedure will possibly lead to an additional source term on 
the right-hand side that accounts for the coupling. To achieve this, we construct a linear 
operator £ which convert two-way fields F to one-way field constituents W, by 

F = CW . (2.22) 

We require that £ when introduced into (2.10) gives, 

C{Id3+V)W = -{d3C)W + N (2.23) 

so as to make V, defined by 

AjC = CV , (2.24) 

a block diagonal matrix of operators. We call £ the composition operator, and W the 
wave matrix. The elements of the wave matrix represent locally the down- and up-going 
constituents. The expression in parentheses on the left-hand side of (2.23) represents the 
two so-called one-way wave operators. The first term on the right-hand side of (2.23) rep- 
resents the scattering due to variations of the medium properties in the vertical direction. 
The scattering due to variations of the medium properties in the horizontal directions is 
contained in V and, implicitly, in £ also. 

To investigate whether solutions {£, V} of (2.24) exist, we introduce the column ma- 
trices, or generalized eigenvectors, £=*=, according to 

£ = ( £+ £- ) . (2.25) 

Upon writing the block diagonal elements of V (generalized eigenvalues) as 

v^(t s- )^ P^"') 

Eqn. (2.24) decomposes into the two systems of equations 

AjC^ = jC^S^ , (2.27) 

where are 2x2 matrices. The central problem that we consider in the present paper is 
to show that there exists an operator pair, {£, V}, such that the above operator equation, 
(2.24), is satisfied. Since the operators are unbounded we need to modify (2.24) and (2.27) 
with respect to the domain of the respective operator. 
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Note that the upper 2x2 matrix of the operators C combines the transverse electric 
field strength, and the lower, the magnetic field strength, whereas the elements of W may 
be physically 'non-observable'. 

We now focus on the fundamental question: Does there exist a composition operator 
£ that decomposes A in the above sense? To begin to show that there exists such a 
decomposition of A, we derive properties of the resolvent set of A which enable us to 
define a certain operator which commute with A. 

3 Properties of the As,x operator 

In this section wc show that the directional decomposition of the electromagnetic field is 
closely related to the spectral properties of the operator A. The definition of the splitting 
matrix requires that there exists a region around the imaginary axis which is free from 
the spectrum. We therefore state the definition of the spectra explicitly. Consider first 
the operator As,x defined on {L^, (•, •)o} with domain H^: 

As,x=A-XI : ^L^ . (3.1) 

Following [37, §5, p.253], [29, §6.5, p.412] and [41, §VIII.l, p.209], we define the spectrmn 
of A, for fixed s as follows: if the scalar A e C is such that the range of ^^.a is dense in 
{L^, (•, •)o} and As.x has a bounded inverse, A is in the resolvent set, P(-4), of A, and we 
denote this inverse by A~\ and call it the resolvent (at s, A) of A. All complex numbers 
not in the resolvent set form a set T,{A) called the spectrum of A. 
To simplify some of the upcoming calculations we use the notations 

s = Sre"^ = Sr cos CT + iSr siu CT and A = Ar + iAi . (3.2) 

The following proposition gives as a corollary that there exists a strip that belongs to the 

resolvent set of A. 

Proposition 1. Let As,x = ^ — A, be defined through (2.14) and (3.1). Let 
Q := {{s. A} e : Re{s} > and (Re{A})2 < (Re {sD^ei/ii} , 

where ei and fii are defined in ^2.1. Then for {s, A} e Q, As,x 

1. is hounded from below; 

2. is one-to-one; 

3. has dense range; 
4-. is closable; 

5. has an inverse: 




where 

£e,x ■■= Ai2 - [An - X)A^^\A22 - A) . 
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Remeirk 1.1. The Hilbert identity or the (first) resolvent equation, 

-4.-i-A-v = (A-A')A-U-v> 

for {s, A}, {s, A'} e Q. // the operator is dosed and if {s, A} S Q then the resolvent is well 
defined and it is an analytic function of X (cf. [24, %in.6.1, pp. 172- 174.], [41> 
pp. 211-212] and [1, p. 84 %3.7.5]). 

Remark 1.2. The underlying requirement of self-adjoint material parameters can he re- 
placed by positivity of the real part of the eigenvalues of the tiuo matrices Sji, se in the 
case of up/down symmetric materials, i.e. when fx^a = l^az = and e^a = CaS = 0, for 
a = 1,2. This follows directly from the proof of part 1. It is not clear that this extension is 
valid for the general anisotropic case, we do not pursue this since our proof of Proposition 
2 makes use of the self adjoint property of e, fi. 

From Proposition 1 it directly follows that: 

Corollary 1.1. For any fixed s G C such that Re {s} > the resolvent set of the electro- 
magnetic system 's matrix contains the strip of all X gC such that 

{Re{X}f <{Re{s}feifii . (3.3) 
That is, the strip belongs to the resolvent set, P(.4). 

Remark 1.3. Let A* be the adjoint of A with respect to the standard inner product in 
L^, denote the spectrum of A by T,{A), and recall the relation T,{A) = T,{A*). Let 

J={'j i), (3.4) 

where I is 2 x 2 unit matrices. From (A. 13) we note that A satisfy the identity 

A{s) - A = J'\A{s)yj* - A = J{-{Ais)y - X)J* , (3.5) 
and consequently when X G 'E{A{s)) then —A G E(yl(s)). 

Proof of Corollary (1.1). Given Re{s} > 0, and the strip (3.3) we find such (s,A) G 
Q, and hence A satisfy the properties 1-3 of Proposition 1. These properties are the 
conditions needed for a point A to be in the resolvent set. □ 

3.1 Properties of a quadratic form 

To prove Proposition 1 we introduce some properties on an auxiliary quantity, a quadratic 
form, Cs^a[^], defined in Lemma 1.1. Using the notation introduced in Proposition 1 and 

in the definition of the norm in (2.21), wc have: 

Lemma 1.1. Let {s, A} G Q and define the quadratic form 

Cs,x[F] ■■= [ se„0E^E0 + s-Via' l(V x E)3f 

+ svc^pHo^Hp + s-h^i |(V X H)3f - 2ARe {E1II2 - E2H1] + 

2i/i33^Im {(V X i;)3M3ai5'a} + 2ie33ilm {ea3^a(V X H)z] d?x' . 

where F = {E, H)^ ■ Then Cg^x is well defined for F G and 

\Cs,y.H*[F]\>C,{s,X)\\F\\l, 

where Cq{s, A) > for {s. A} G Q. 
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Proof. That Cs,\ is well defined for F G is clear as the quadratic form contain at most 
one derivative for each field component in each term. To see that Cs,\ is bounded from 
below we take the real part of the integrand and using the notation introduced in (3.2) 
we obtain 

{SrCOSa)So:0EaEj3 + {s~^ COSa)lJ,^^ |(V X E)3\^ + [s^ COS u)VolI3H ail p 

+ {s-^ cosa)e33^ |(V x H)z\^ - 2ARRe - E2H1} , (3.6) 

due to that s and u are self-adjoint. Since {s, A} G Q, it follows that costr > and we 
have 

(3.6) > cos a {e^pE^Ep + UccuH^H,}) - 2\^Ke {E1H2 - E2H1} . (3.7) 
For all ry > we have the inequality 

2ARRe \E1H2 - E2H1} = 2ARRe \ E^H \ < r?ei|^PsrCosc7 + — — ^ 1^1^ , 

^ ^ J ryeiSrCOSCT 

which implies 

(3.7) > eiSrCOsa(l-r])\E\^ + ( fiis^ cos a - ) . (3.8) 

\ rjeiSrCOsa J 

Thus we require that t] G (0, 1). The largest |AR|-strip is obtained in the limit — > 1, thus 

A^ < (Re{s})2eiAi . (3.9) 

Hence for given fixed s such that Re{s} > and for a fixed A that fulfils (3.9), there 
exists an optimal 77 such that the bound from below in (3.8) is maximal. Thus the best 
choice of bound from below with the given estimates is 

Cq(s,X)= max min fRe{s}ei(l - r?), f/iiRe{s} ^ ^ > . (3.10) 

o<r,<i V V r?eiRe{s}y/ 

We note that the argument above is a continuous function of r] on [0, 1], and the maximal- 
value is attained (no need for sup) in the interior of the interval. Furthermore, the optimal 
value T] for which the Co(s, A) is obtained as a solution of a second order equation in 77, 
but here it suffices to know that it exists and is positive. Upon integration we find that 

\Cs,x[F]\ > Co{s,X){\\E\\l + \\H\\l) = Co{s,X)\\F\\l . □ 
3.2 Proof of Proposition 1, part 1 

To start the proof that the operator As,\ is bounded from below, we employ Schwartz' 
inequality 

llA.Ai^llolli^llo = \\KAs,xFMF\\o > \{F,KAs,xF)o\ , (3.11) 
with F G D(^) = and where 

^=(? ' ^'-''^ 

and / is a 2 X 2-matrix. K is a, unitary for || • ||o. On the block element level where 
F = {E,H) {cf. (2.12) we have 

Fg{K{A - X)F)q = E*A2iE + H*Ai2H - 2ARe + E*A22H + H*AiiE , 
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where the repeated index q indicates summation over the four components and E* = 
{Ey . The first term becomes after integration over x' and integration by parts, 

/ E*A2iEA^x' = f s£c«/3^ai^/3 + s"V33M(Vxi;)3|^ dV (3.13) 

and the second becomes 

/ H*Ar2Hd?x' = [ su„0Hc,H0 + s-'^e^^\{V X H)3f d^x' . (3.14) 
The fourth term becomes after simplification and integration by parts 

/ E*^22^dV = / e^3le„3^„(Vxif)3-(Vx^)3A*33V3a^fadV , (3.15) 

and the last term 

/ H*AuEd''x' = [ (VxE)3Mi3Va3^a-ei3'e3a^a(Vx5-)3dV . (3.16) 

Since e and are self-adjoint matrices the two terms in (3.15) and (3.16) combine to 

/ 2iIm|^M22^| dV = / 2in^^Im{{V X E)3fMa3H„} 
Jr2 J 

+ 2ie^^Hm {e„3^a(V x H^} d^x' . (3.17) 

Thus, 

\iF,KAF)o\ = \C,,x[F]\, (3.18) 
for F e D(^s_a) C By Lemma 1.1 and (3.11) we obtain 

\\AF\\o>Cois,X)\\F\\o , (3.19) 

where Co(s, A) is defined in the lemma, and Co(s, A) > for {s, A} G Q. □ 

3.3 Proof of Proposition 1, part 2 

The inequality 

\\As,xF\\o>Cois,X)\\F\\o, (3.20) 

with Co > from part 1, implies that the null space only contains the zero element. 
By [29, p. 171, theorem 4.4.1] an operator with trivial null space is one-to-one (injective). 
Hence, the operator As,\ is one-to-one for {s. A} e Q. □ 

3.4 Proof of Proposition 1, part 3 

Let A* ^ denote the adjoint of As,x on L^. To show that the operator has dense range it 
is sufficient to show that the kernel of Al x is trivial. That is, 

AlxG = , G € D(^*) ^ G = , (3.21) 

thus if Al X is bounded from below then the desired result follows directly, cf. Proposi- 
tion 1, part 2. 
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The adjoint of As,\ with respect to the inner product (•, •)o is 

11 ~ -^21 _ 

-4*2 -^22 ~ 



AU = { ./is, ) , (3.22) 



where 



j^* ^ _ ( ^2^23 ^2^J.13 \ -1 _ ^-1 / eai^i -631192 
" V ^l'"23 dip,i3 J 33 -63291 £32(92 

A* =s( ''J'^ ^^"^ \ - s-^ ( ^^^^^^^ -die-3id2 
V ^^21 J^ii y ' V -d2e.yldi i92e33^i92 

A* =s( ^ - S"^ f ^2^33^92 92Mi3^5l 

.* _ _ -1 ( M3292 Ai329i \ 9iei3 -91623 , -1 

•^22— A*33 r,„.^„ r,„,Pl, -Fl„a,„ Fl„-S„„ ^33 ' 



^31^2 J \ -^2613 ^2623 

where we have used that e, are self adjoint, and hence that their diagonals are real- 
valued. The domain of the adjoint is the set 

D(^^) = {GgL2: 

3H e L2 such that (if, F)q = (G, As,xF)o VF € D(A,a) } • (3.23) 

To show that A* ^ is bounded from below we will use the same method as in Proposition 1, 
part 1. First we need a small enough set that contains the domain; from the form of A* )^ 
it is clear that if G = A* e L^, then at least F e H^. Thus the domain is contained 
in H^, that is 

D(^:, J c Hi . (3.24) 

To obtain the quadratic form needed to use Lemma 1.1, we once again use the Schwartz 
estimate. Let G G D(^* and F := KG, where the matrix K was introduced in (3.12) 
and is unitary on (•, •)o and has the properties K = K* = and ||ii'G||o = ||G||o. 

Similarly to (3.11) we find 

IKaG||o||G||o = IKaGIIoII^^IIo > \{F,Al.G)^\ = \{F,AUKFU = 



\{KAs,xF,F)o\ = \Cs,x[F]\, (3.25) 
and with F — {E, H}. Let {s, A} E Q, from the properties of Cg.x, we find 



\{F,Al^KF)o\ = \Cs,x[F]\ > Co{s,X)\\FU = Co{s,X)\\GU , (3.26) 

where Cs,x is defined for F e H^. From (3.25) and (3.24) it follows that A^ ^ is bounded 
from below for all G e D(^* _^) c H^. For the bound from below of A^ ^ it directly 
follows that A* ^ has trivial kernel and thus As,x has dense range for the condition of 
{s,A}eQ. ' ' □ 

3.5 Proof of Proposition 1, part 4 

In Section 4 below we consider a Dunford- Taylor integral over the resolvent and the 
analysis simplifies if A is closed. Prom the form of A* ^ we note that C°°(M^,C^) C 
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DiA* ^) and hence it is densely defined for {s,A} € Q and thus, by [24, p. 168, §111.5.5], 
the operator As,x is closable. The closure is denoted by {^s^aIci- 
The range for the closed operator is still dense in {L^, (•, •)o} since 

R(^,,a) CR(Mm}J CL^ . (3.27) 

for {s. A} G Q. 

To show that the closed operator is bounded from below, we rely on Corollary VI. 1.19 

[24]. This corollary applies to sesquilinear forms in Hilbert spaces, but due to example 1.23 
and example 1.3 in [24] we draw the conclusion that we can construct the sesquilinear 
form {As,\F,As,\G)o and that it is only closable when As,\ is closable. Thus, by the 
above mentioned corollary, we obtain that the closed form is bounded from below with 
the same constant and thus, the closed operator is bounded from below for {s. A} G Q. □ 



3.6 Proof of Proposition 1, part 5 

Given (s, A) S Q we have shown in part 1-3 that As,\ is one-to-one, has dense range 
and it is bounded from below. Consequently we know that the inverse exists and is 
unique. The operator As,\ can be explicitly inverted in terms of the inverse of two 2x2 
matrix operators through a quasi-diagonalization. Once again introduce the matrix K 
{cf. (3.12)), with K-^ = K. We find 

A _ / ^21 A22 — A 

V -^11-A A^2 
and from the form of A2i, (2.14) we find 



I (^,.421^)0 1 



seocpEc,Ei3 + s-Via'KV x ^)3|^dV 



>ReMei||^||^ 



(3.28) 



hence the inverse exists and is bounded for Re {s} > 0. Thus .421^ is well defined. The 
operator KAs,\ is then diagonalized as follows 



where 



and 



while 



with 



T -( ^ 

V {A^^-X)A:,^ 1 



To 



2;s,A 



V, 



1 -A^l{A22-\) 
1 



A21 
£s,x 



(3.29) 
(3.30) 
(3.31) 

(3.32) 
(3.33) 



£s,x = A12 - {An - X)A2i{A22 - A) . 

The characteristic operator, £s,\, a matrix extension of the of the 'transverse Helmholtz' 
operator [14]. For each fixed A, the operators '7i;s,A)^;s,a have the inverses: 



1 
-{Au-X)A^i' 1 



(3.34) 
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and 

T-, = ( J A-i^(^--A)^ , (3.35) 

respectively. From the quasi-diagonalization (3.29) we obtain an explicit expression for 
A~\ in terms of f ~^ . That f ~^ is well defined follows by considering 

\{E,£,,xE)o\ = \{E,Ai2E)o + iiA22-iXi)E,A2i\A22-iXi)E)o-XliE,A2iE)o\ (3.36) 

where we have used Ali = —A22- The bounds from below A12 > Re{s}/ti and A21 > 
Re {s} ei, gives 

\{E,£s,xE)o\ > R^^llReisD'eiAi " Xl)\\E\\l , (3.37) 

and since (s, A) G Q the bound from below is positive and £s,X is invertible. Starting from 
(3.29) and inverting term by term, gives 




(3.38) 



□ 

4 The Splitting Matrix 

We proceed with the decomposition of the electromagnetic system's matrix. As the spec- 
trum is absent from the strip (see Corollary 1.1), we define a certain commuting operator 
through a resolvent integral, this operator will satisfy a number of properties, and will 
be called the splitting matrix. Wc note that if an operator has a spectral resolution, or 
even a part of the spectrum which is bounded, then one can define a projector with help 
of a Cauchy type integral, also called Dunford's integral, over the resolvent with inte- 
gration path around the bounded spectral region, sec [24, III. 6. 4], and also [41, 32]. For 
the electromagnetic system's matrix such information about the spectrum is not known, 
we do know however that the spectrum is separated into two parts by a strip around 
the imaginary axis. The idea here is to accomplish a decomposition by introducing an 
operator defined by an integral over the resolvent of A, as to try to split the two parts we 
know exist, similarly to the case of bounded spectral regions. We use the Dunford- Taylor 
integral applied to a closed, imboimdcd, operator as in [4], cf. [24] for accretive operators. 
This theory is given only for closed paths, or absolutely bounded integrals, hence the 
extension needed here to non-closed paths require that we prove that the operator is well 
defined. 

In this section we prove a number of properties of the splitting matrix, among them 
that it is well defined as a pseudodifferential operator with a parameter, that it is an 
involution, and that it commutes with the electromagnetic system's matrix. Once the 
splitting matrix is shown to be well defined we derive its generalized eigenvalues and 
eigenvectors; the generalized eigenvectors are the key components for the decomposition 
detailed in the next section. 
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4.1 Definition of the Splitting Matrix 

Given a fixed positive constant S'r > 0, let 



Qi = {{s,A} e : Re{s} > 5r, |args| < tt/2 and |Re{A}| < SnVhjh} ■ (4.1) 

From Proposition 1 we note that Qi C Q, hence by Corollary 1.1 the strip |Re{A} | < 
Sny/eipa belongs to the resolvent set of the electromagnetic system's matrix. Thus we 
can consider the operator defined through 

B= lim \ I dA {{Aj^^-IXy'^ . (4.2) 



n— »oo TTl 



cl 

AeK„ 



The spatial and time-Laplace dependence is present but not explicit in the notation. The 
integration path is: 

K„ = {A e C : Re {A} = t/2 and |Im {A} \<n} , (4.3) 

where 

T = S'r^/IIAi , (4.4) 

and hence the integral path is in the resolvent set. Some of the considerations that follow 
become simpler if we consider the operators restriction to C°°(M^,C^), and similarly for 
any operator with the notation . With the above introduction, we have the following 
proposition. 

Proposition 2. Let B be defined as in (4.2) with {s,K„} e Qi then B 

1. is a pseudodifferential operator with parameters of order 0; 

2. has a restriction B\^, which maps into H*~^; 

3. 'commutes' with A in the sense that on the set H'^ we have 

B\-^A = AB\^ ; 

4- has a restriction, B, that is an involution; 

o o 

5. has a restriction B that has a generalized eigenvector ; unique up to a normaliza- 
tion, and with a corresponding scalar 'eigenvalue' 7 = ±1, satisfying the equation 

BC^ = 7£± . 

o _|_ 

The explicit form of is 

^±^({±i + Bi,)mA 

where is a normalization in the form of invertihle 2x2 operator matrices; 

o o 

6. has a restriction B that is one-to-one on a core and thus its element B21 is invertihle 
on its range; 
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7. has 'generalized eigenvectors', i.e., that is the extension of exists. With proper 
choice of normalization, £±1^ : H«(IR^C2) ^ H«-i. 

o o — 1 o ~ 

Remark 2.1. With the choice of Af as B21 , we identify C i as a mapping between H 
and E, i.e., an impedance mapping. The corresponding map in linear acoustic is a map 
between the pressure and the vertical particle velocity cf. [22]. Both these mappings are 
the acoustic, and electromagnetic respective equivalent maps to a Dirichlet-to-Neumann 

map for the wave equation [11]. 

Remark 2.2. To ensure that the above defined splitting matrix is non-trivial we have to 
exclude two cases: that the integral (4.2) collapses to the identity or to the zero operator. 

Whether this happens depends on the non-triviality of the spectrum. To ensure that this 
is not always the case we consider the homogeneous-isotropic case, see Appendix B where 
we obtain the explicit form of B. We note that if the medium in the neighborhood of a 
point is isotropic, then microlocally, at that point, the operator B reduces to the isotropic 
case, which is clearly different from the identity and the zero operator, and hence B can 
not be the unity or the zero operator for such media. 

4.2 Proof of Proposition 2, part 1 

The operator B is defined through an improper integral over the resolvent. To prove 
that B is well defined as a pseudodifferential operator with parameter we consider first 
the paramctrix of As,\ and then integrate each term of the asymptotic expansions with 
respect to A and prove that this step is well defined. Hence we obtain an asymptotic 
expansion for the symbol of B, via the usual calculus of pseudodifferential operators we 
thus construct a well defined operator B. 

4.2.1 PseudodifFerential preliminaries 

The calculus of pscudodifFcrontial operators can bo introduced by means of a Fourier 
transform, thus defining signs and symbols. For simplicity we use standard notation for 
the symbols and their compositions. TliroTighout this paper we use the left symbol (in 
the notation of [32]). The Fourier transform, in the plane with respect to the first two 
variables x' = {xi,X2} has an inverse given by 

E{x,s) = -^[ e^^'-^'{J^E)ie,X3;s) (4.5) 

for the complex field E e L2(R^C2). Here C ■ x' = £,1X1 + £,2X2. To obtain the left 
symbol of A21 we let it act upon (4.5) and obtain that the integrand expression in front 
of e'«'-^'(jr^)(^',a;3,s) is 




This is the left symbol of ^21 • To find the appropriate behavior of the symbols we have 
to consider symbols with parameters. We consider s to be a parameter of the same order 
as We hence find the principal symbol of 021 to be 



B. L. G. Jonsson 



2008-09-04; 11:22 



Wave splitting of Maxwell's equations. 



17 



and a2i;i is homogeneous of order one in (^',s). 



4.2.2 EUipticity of As,x 

For matrix valued operators it is the determinant of the symbol that controls the regularity 
and existence of its parametrix. Wc require that the coefficients to {£,', s, A) arc arbitrarily 
smooth for each term in As,\, and thus we can use the criteria in Definition 5.1 together 
with Proposition 5.1' of [32, pp. 38, 39] to define ellipticity of As,\- To construct the 
principal symbol of the operator A, a-i{x,^';s), we proceed as above and obtain for the 
remaining elements 

ailil — IM33 t t I + 1*^33 t f > K^-') 



-1 / M32?2 Ai31^2 \ I • -1 / £13^1 ~ei3^2 



• — i I M'S^^.Z f^JiS2 \ I • — i / '-iJSi >-lJS^ \ I A c\\ 

"^^^^ = \ m.i. M3l6 ) + \ -£236 £236 ) ■ 

Let 

OL-x ■■= a.\ - A/. (4.10) 

Then a.-^ have homogeneity degree 1 in (6,s,A). The remaining part of the symbol of 
As,\ has a lower degree of homogeneity in (6, s, A). To ensure the ellipticity of an operator 
in the parameters (61 s, A) the following estimate is needed 

Ci(ie? + |s|' + |A|2)2< |deta;i| <C2(|6|' + |s|' + |A|2)2 (4.11) 

for some R such that |6P + |s|^ + |A|^ > i?^ and with proper restrictions on the parameters 
{s,A}. The properly supported requirement for ellipticity follows from the fact that A 
is a classical pseudodijferential operator with smooth coefficients [32]. The upper limit 
of (4.11) follows directly from the fact that deta;i is a polynomial, homogeneous of 
order four in (6,s,A) (sec Appendix C), together with the fact that we can dominate 
this polynomial by (|6P + + ^^"^ ^ constant, for some constant Rq such that 

I6P + 1^1^ + l-^P > -Rq- To prove the lower limit we need a more subtle method. 



The lower limit of (4.11): This is a multi-step proof, and a complication arises since 
the region in s, A is not conical. First we prove that the determinant is non-zero on a 
surface (sec Figure 4.1), then we use a scaling argument to extend this to a bound from 
below of the form (4.11) for a conical region with the surface in figure 4.1 as 'bottom 
surface'. In the last step we extend the obtained result to the non-conical domain, so as 
to include Ar. 



Non-zero determinant of CC;! Let Ar = 0, to explicitly show that the determinant 

of a-i is non-zero is difficult due to the large number of terms that it contains (see 
Appendix C). Our scaling argument needs only that the determinant is non-zero on a 
surface, here part of a sphere, see figure 4.1. Thus let |s|^ -|- |^|^ -|- jAij^ = some 
positive constant R^. Wc consider two cases; s = and s ^ 0. For the first case with 
s = 0, we find from Appendix C that 

deta;i|,^o = 

(Ai - Aie33^(ea3 + e3a)6 + £33^£a/366) (''^I " '^I/^SsHMtS + fJ-S-y)^! + ^33^7566) • 

(4.12) 
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Figure 4.1: A schematic picture of the surface + |sp + |Aip = R^, |args| < 7r/2. 
In particular note that the condition | args| < 7r/2 shrinks the surface away from a half 
sphere. 

Using the restriction that /i and e are self adjoint, together with the estimate 

2Ai^ae33^Re {es^} < rj^'^Xi + Vi^ai^asea/i^a^i} (4.13) 
for rji > and similarly for 2Ai/i33^Re {/U^s} ^-y gives 

I deta;i|,=o > ((1 - %-^)|Aip + Ms^s'^S^C.?/?) ((1 " + e^s'^^^l^^si.) 

> Ca (|Al|2 + le?)' = Cai?^L=o,A«=0 ' (^-l^) 

where i?2|^^^_^^^^ = |Ai|2 + |^f and 

^a/j' = ~ ^lecseia^es/J and u^^^^ = iiap - V^l^aSjJ'si 1^3/3 ■ (4.15) 
To find an explicit expression for Ca, let us choose r/i > 1 so that 

< 1 - ryfi = inf e^^ inf e^'kj^S • (4-16) 

x' K 1=1 

The right-hand side of this equation gives the minimum of the lower eigenvalue of the 
matrix e['^g '' normalized with €33. There exists an r^i > 1 that fulfils this equation since 
ei > 0, hence the left-hand side of the above expression is positive. The same way we find 
an > 1 such that 

< 1 - 7,2^1 = inf M3-3' inf ^.(7)^76 • (4.17) 

x' |{'|=1 ' 

The constant Ca becomes 

= (1 - %-')(! - % ') • (4.18) 
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For the case where s 7^ we use Schwartz' inequaUty on an inner product. Thus we 
introduce the 'matrix' norm 

4 

\F\s=y2\Fi\'' and |a;i|s = sup \oc.,iF\, , (4.19) 

1^1 l^ls = l 

with a corresponding inner product defined analogously and denoted by (•, •)s- Both the 
norm and inner product depend on (x, s, Ai, 0). Consider the normal 4x4 matrix 
a*ia-i, for Ar = 0, that have eigenvalues ki,...,K4, each with a variable dependence 
(x, s, Ai, 0). From the definition of eigenvalues it follows that < Kj S R, where 
i = 1, . . . ,4, and we use the convention > • • • > ki. From the relation 

Ki < K4K3K2K1 = det Oi*ia-i\ = |deta;i|^^^g , (4.20) 
we find that it is enough to prove that Ki ^ 0. Schwartz' inequality (c/. (3.11)) gives 

|a;li^|s|i=^|s>|(i^,i^a;li^)s| . (4.21) 

Thus if for Ar, = and s ^ we can obtain an estimate of the form 

\{F,Kcx.,iFU>Ct\F\', , (4.22) 
where Cb > 0, then from (4.19) and (4.22) it follows that 

^-1/2 = < (O-i i^gjjpg Ki>C^ . (4.23) 

By (4.20) we obtain, 

\deta,i\,^^^>Ct , (4.24) 
under some restrictions on s, to be derived. Now with the explicit form of oc-i we obtain 

{F, Ka;lF)s\^^^Q = S {sa^EaEff + Ua^HaH^) 

+ s"' (m3"3 I(C X E)sf + |(C X H)3f) - 2iAiRe {EJI2 - E2H1} 

+ 2i (Mia^Im {(e X E)3i^a3Ha} + e^^Im {e^sEai^ x H)^}) , (4.25) 

where we have used the notation of Proof of Proposition 1, part 1. Since s 7^ we take 
the real part and obtain 

\{F,Kcx.,iFl\x^=o > Re{s} {ei\E\^ + ili\Hf) > Cb\F\l , (4.26) 

if Re {s} > and here 

Cfo = |s| min{ei,/ti}cos(T (4.27) 

is positive if |c7| = | args| < ■jt/2 and s ^0. By the above argument, (4.14) and (4.24) the 
determinant is non-zero on the surface 

Rl = l^'\ + \sf + \X,W |(7| <7r/2 (4.28) 

if Re ^ 0. Hence there exists a lower constant Cg = min(C(,, CaR^) such that 

|deta;i|^^^o>Ce>0 (4.29) 

on this surface. 
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Figure 4.2: A schematic picture of the conical region in s, Ai. Note that the outer 
surface can have an arbitrary radius larger than the inner radius Rg. 



A scaling argument: To extend the result 

det a-iix, £,'; s^, a, Ai, Ar)|^^^o 



> Ce > 



(4.30) 



for l^'l + |sp + |Aip = i?g and |(t| < 7r/2 to a proper bound from below, wc use a scaling 
argument. The homogeneity of deta;i allow us to scale ^',s, Ai, to an arbitrary radius 
greater then R^, and 

det e'; s,, a, Ai, 0) = i^'^'d^T + \s? + | Ai^^ det cx,i{x, f; Sr, cr, Aj, 0) , (4.31) 

where the ~ variables are normalized to lie on the surface + |sp + |Aip = R^. Thus 
from (4.30) we have obtained 



|dcta;i(x,r;Sr,a,Ai,0)| > C,R-W$,'\^ + \s\^+\\i\''f , 
in the conical domain, 

^' e , Ai e R , s e C and | args| < 7r/2 , 

see Figure 4.2. 



(4.32) 



(4.33) 



The case Ar ^ 0: To extend the argument above to include the case Ar, ^ we 
impose the condition Re{s} > 5r — ^SnV^ifi-i- Let Ar < r and Rg > Sr. For large 
enough Re, where |^'|^ + |sp + |Aip > R^, the worst case for the bound from below of the 
determinant is 

I det a;i| > CeR^^ (ler + \s\^ + \Xi\^) - CrT + l^l^ + lAi^ (4.34) 
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where Cr is chosen to include the sum of the maximal material parameters in front of Ar,. 
Using the Holder and the Jensen inequalities [10, p. 28, Theorem 19] gives 

I det a;i| > CeR-' (ler + \s\' + lAiH - CrT + \s\' + \Xi\') 

> CeR-" (iri' + Isl" + |Ai|4 + |Ar|4) - CrT + \s\' + \Xl\' + M) 

> CeR-^4-^\z\^ - CrT\zf = A-^C^R'^lzf {\z\ - 4%^t) , (4.35) 
Cr = max{Cr, CeR~^} (4.36) 



where 
and 



\z\' = \e\' + \s\' + \X\' . (4.37) 
Comparing with (4.11) we obtain the condition 

Cl\z\'^ < 4-^CeR-^\z\^ - Cr\z\^T . (4.38) 

Thus for some given, arbitrary Re > 0, there exists a large enough R such that | det cx-i\ > 
Ci\z\'^, for \z\'^ > i?^ when (Ci,i?) satisfy the following constraints: 

< Ci < 4:-^CeR~^ (4.39) 

and 

\zf >R^> ma^iRlRl + r^ . (4.40) 

Hence for ^' £ and {s, A} G Qi, the operator is elliptic. 

Note that the quadratic form argument of Lemma 1.1 can be applied to the symbol 
cx-i to yield a positive lower bound. Consequently (4.23) and (4.20) imply 

|deta;i|^0, (4.41) 

for {s, A} e Q and \s\^ + + |A|^ ^ 0. However the desired increase in |s|^ + + |Ap 
does not follow directly, since the domain is not conical. Observe also that (4.41) is true 
in the region {s, A} S Qi. This result will be used in the end of the proof of part 1. 



4.2.3 The parametrix of As,\ 

We have above shown that ^^.a is elliptic in pscudodiffcrcntial sense, hence the corre- 
sponding parametric is well defined. In the subsequent analysis we are interesting only 
in the principal part. Prom Proposition 1, part 5, we know that the inverse can be effi- 
ciently expressed in terms of £~j^ for s, A e Q and different combinations of 2 x 2 matrices. 
Therefore we introduce the notation of ^ on 2 x 2 matrices defined by 

a2i;i := ( ^''''■''^f -Ki;i)i2\ (442) 

V^('^21;l)21 (a21;l)ll / 

where is the (ij)-element of the 2x2 matrix. From the definition it follows directly 
that 

a21;l = a21;l , (flll;! 021;! ) ' = a21;iaii;i , (4.43) 

(aii;i + a22;i)^ = aii;i + a22;i , (4.44) 
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and 

a2i:ia2i;i = dcta2i;i/ . (4.45) 
The principal symbol of the characteristic operator, £s,\, is 

e;i = ai2;i - (aii;i - /A)d2i;i(a22;i - /A)(det a2i;i)"^ , (4.46) 

and using (4.43)-(4.44) we find 

e;i = ai2;i - (a22;i - /A)a2i;i(aii;i - /A)(det a2i;i)~^ , (4.47) 

thus 

(dete;i)/ = e;ie;i . (4.48) 
Prom writing out all terms we find that 

deta;i = (dete;i)(deta2i;i) , (4.49) 

which is a polynomial homogeneous of order 4 in s, A. That | det e-i \ ^0 follows directly 
from that | det a-i \ ^ (sec (4.41)) together with the observation that | dct a^i-.i \ ^ and 
that I deta2i;i| is bounded above and below by constants times |^'P + |sp. Hence it follows 
that the parametrix of £s,\ is well defined. This is to be expected since the inverse of 
£s,\ was shown to be well defined in Proposition 1. part 5. With the above consideration 
we find that the components of the principal symbol of the resolvent, T-^-i ;= (y~^ = 
(a — A/)~^, are (c/. Proposition 1 part 5) 

(r;_i)ii := -(deta;i)"^d2i;i(a22;i - lX)e.i , 

(r;_i)i2 := (deta2i;i)"^a2i;i(/+ (deta;i)"^(a22;i - /A)e;i(aii;i - /A)d2i;i) , 
(r;-i)2i := (deta;i)-ie;i(deta2i;i) , (4.50) 
(r;-i)22 := -(deta;i)~^e;i(aii;i - /A)a2i;i . 

Upon integration of the parametrix with respect to A, the element (r;_i)i2 has an unsuit- 
able form, therefore we use the identity e-ie-i = 7dete;i and rewrite (r;_i)i2 into 

(deta;i)(r;_i)i2 = d2i;i(e;ie;i + (deta2i;i)~^(a22;i - /A)e;i(aii;i - /A)d2i;i) 

= d2i;i(e;iai2;i + (deta2i;i)~^[a22;i,e;i(aii;i - /A)d2i;i]) , (4.51) 

where [•,•] is the standard commutator [A,B] = AB — BA. By inserting the explicit 
expression of e;i in the commutator together with the two relations (dii;i — /A)(aii;i — 
/A) = /det(aii:i — /A) and d2i;ia'2i:i = dcta2i;i/ we find that (4.51) reduce to 

(deta;i)(r;_i)i2 = d2i;i(e;iai2;i + (deta2i;i)"^[a22;i, di2;i(aii;i - /A)d2i;i]) . (4.52) 

An alternative form of (4.52) is obtain by inserting the explicit form of e-i and simpliiying: 

(detQ;;i)(r;_i)i2 = d2i;i(/det ai2;i - (det a2i;i)~^(d22;i - -fA)a2i;i(dii;i - /A)ai2;i 

+ (deta2i;i)~^[a22;i,di2;i(aii;i - 7A)d2i;i]) . (4.53) 

With the above expression we have obtained the principal part of the symbol of A~\ 
and each term in the matrix has the form 

A" 
deto;;! 
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times a constant homogeneous in ^, s of order 3 — n, here n = 0, . . . , 3. For the lower order 
terms we use the recursive construction formula of e.g., [32, pp. 44, 45 §1.5.5] to deduct 
that symbols of lower orders have the A dependence 

A"' 

(deta;!)™ ' 

where 4m — n' > 1. We have above deduced the A dependence for all terms in the symbolic 
expansion of the symbol of ^^A' furthermore, we have the principal part explicitly. 

4.2.4 ,8 is a pseudodifferential with parameter of order 

Given the parametrix of the resolvent, we integrate each term of the asymptotic series with 
respect to A. To validate this procedure we show below that each of the terms is finite and 
that the integration docs not rearrange the terms with respect to order, i.e., the principal 
term remains the principal term. We also show that B is an operator corresponding to a 
symbol that is homogeneous of order in s). 

As shown in the previous section, each element of the resolvent has the form 

A" 

(deta;i)™ ' 

with Am — n > 1. The principal symbol has m = 1 and n = 0,1,2,3, all other terms 
have homogeneity or lower. In the evaluation of the integral we distinguish between two 
different cases, the principal valued integral corresponding to n = 3, m = 1 and the other 
cases. We observe that due to the homogeneity of the parametrix terms the case n = 3, 
m = 1 is the only principal integral. 

The case m = 1, n. = 'd gives a finite result, which we find by evaluation the integral 
over the integrand (4.55). In order to do this claim that we can use the representation 

det a;i = (A - Af )(A - Aj)(A - Ar)(A - Aj ) , (4.55) 

where the eigenvalues, i.e., the roots of the fourth order polynomial deta;i, are denoted 

by Af , A^ where the +(— ) indicates that they have positive (negative) real part. Indeed, 
to show that two of the eigenvalues of a-^i have positive(negative) real part we consider 
the isotropic case. The isotropic deta;i have the A-roots {of. Appendix B) 

X = ±^/s^^^^°ij}^°TW , (4.56) 

i.e., two double roots on each side of the strip |Re{A} | < Samix' \/ e'^°/i'^°. The lower 
bound on As,\ in Proposition 1, part 1 shows that for each anisotropic material with 
instantaneous response, the area around the imaginary axis is free from eigenvalues. We 
introduce a parameter 7 in the material coefficients by 

e(7) .= giso + _ giso) (4 57) 

and analogously for /i^'^^ The lower bound on a;i, that is obtained by applying Lemma 1.1 
and Proposition 1, part 1-3, to ol;\, and apparent in (4.41), ensures us that for all 7 G 

[0, 1] the I det af^^ | is bounded from below and that there are no eigenvalues on the 
strip around the imaginary axis, and since the eigenvalues of a matrix depend point- wise 

continuous on its coefficient [24, pp. 107-108, §2.5.1], the eigenvalues vary continuously, 
but not discontinuously on each side of the imaginary axis. In the case 7 = there are 



(4.54) 
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two eigenvalues on each side, by counting multiplicity and hence, by the continuity of 
the eigenvalues, this has to be the case for all 7 e [0, 1]. Thus we have shown the claim. 
The symmetry of A discussed in Remark 1.3 can be used to show the same result for the 
special case s = Re{s}. With the representation (4.55) we evaluate the integral 

/n J.3 
dAi , (4.58) 

_„ det ct;! 

by partial fraction decomposition. Assume initially that there are no equal roots then 



deta;i A-A;;^ A-Aj A-Ai A- A 



(4.59) 



2 



where all D depended only on A^2- Each such fraction is integrated over the imaginary 
axis to become 



lim / 



dAi lim ln(A- A^2)k,__„ 

A-Aj^a " (4.60) 

= i7rsgn(Re{AR - Xf ^}) = Titt , 



where the branch cut is along the negative imaginary axis and where |Ar| < |Re {A^j} I- 
Concerning the choice of branch cut, observe that the integral above should be summed 
over each eigenvalue, thus the branch cut of the logarithm has to be chosen such that it 
agrees for all eigenvalues, hence the negative imaginary axis. Thus 

/" A^ 
dAi XTT— = 17^(^1,+ + ^2,+ - Di- - D2,-) (4.61) 
-n 



n— ^00 



det a-i 



for n = 0, . . . , 3. Here 

Di,- 
D 





(A^ 


-A+)(A+-Ar)(A+ 






-(A^)^ 




(Af 


-A2+)(A+-Ar)(A+ 


-A2) 




(Ar)' 




(Af 


-Ar)(A+-Ar)(Ar 


-A2-) 








(Af 


-A2)(A2+-A2-)(Ar 


-A2) 



2,- •— 

and each term is homogeneous of order 0. To show that the sum is bounded from above 
we have to eliminate (A^ — A^) and (A^ — A2 ) from the denominator. We find that 

n A+A+(A+ + A+) - Ar((A+)^ + A+A+ + {X+r) 

""''^ ^ ""''^ - (Af-Ar)(A+-A2-)(A+-Ar)(A+-A2-) ^''''^ 

and 

(A+-Ar)(A+-A2-)(A+-Ar)(A+-A2-) ^ ' 

hence the denominator is bounded away from zero by (2t)*, since |A^| > r = SR\/e\(ii, 
this follows from the Corollary 1.1 that shows that the strip |Ar| < r is free from eigen- 
values. Hence the integral is bounded. 
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The case with equal eigenvalues follows similarly. Assume A2' = Xf and ^ A2 then 

A^ _ L»i.+,2 D2- 
deta;i " (A- A^)2 ^ A- Af ^ A - A^ ^ A - Aj ' 



where 



(Ar-Ar)(Ar-A2-) 

A+A+(A+A+ - ArA2) + (A+ + A2+)(A+ + A+ArA2- - (Af - A2-)A+A+) 



(A+-ArF(A+-A2-)2 
With the integral 

Xi—n 



r 

lim / 

n—*oo I „ 
t/ —n 



dAi — — r-r = lim 



(A - Xt)^ n^oo Xt -X 

and (4.60) we find that 



(4.65) 



At =— n 



\3 

lim / dAi — = i7r(Di,+,i - - D2,-) , (4.66) 

— ' deta;i 



n— >oo 



and hence, it is bounded from above and homogeneous of order 0. The case where X^ — X2 

and X2 7^ A^ is totally analogous. 

For the case with two equal eigenvalues we have 



dcta;i (A-A+)2 A-A+ (A - Ai")2 A - A2 
and from (4.60) and (4.65) we find that 



(4.67) 



lim / 



X^ 

dAi — = i7r(L»+,i - 7^_,i) . (4.68) 

det OL-i 



Hence we need only to find -D+^i and -D_^i, 

(A+)2(A+-3Ar) ^ _ (A-)2(3A+-Ar) 
— — ■ — — — • ^ 



and hence the integral is homogeneous of order and bounded from above since the 
denominator is bounded from below. We have thus shown that the integral (4.58) is well 
defined, and homogeneous of order for all possible combinations of A-roots in det a-\. 

Next we show that the remaining terms have homogeneous degree +1 compared to 
the corresponding term in the polyhomogeneous expansion of r and consequently that the 
integral over A does not rearrange the symbol expansion. Prom the construction of the 
parametrix of As,\ we know that its asymptotic symbol expansion has a A dependence 
of the form (4.54). Using that the determinant deta;i is homogeneous of degree 4 in 
(5', s. A), we find that 

/A" 
(deta;i(a;,f ;s,A)) 
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is homogeneous of degree n — 4m + 1 if 4m — n > 2. Indeed, 



^n— 4m+l 



dA r] 



{deta.i{x,ri^';r]s,r]X)) 
A" 



(4.70) 



(deta.i{x,r]^';r]s,X)^ 



= In,m{x,ri£!;r)s) 



where we have used 1) that the hmits of the integral goes to infinity, 2) that the A-roots 
scale with r], implying that the strip Ar < ryr is free from poles, 3) the scaled integration 
path is equivalent to the integration path of Ar = r/2 since A~\ is analytical in A in 
the resolvent set and 4) that r/ is such that r/Rc {.s} > S-r,. Thus In,m is homogeneous 
of degree n — Am + 1 in s. Each term in the poly homogeneous expansion of r have a 
A-dependence in the form of In,m with a A independent coefficient. It follows that each 
integrated term of the expansion has a homogeneous degree that is one order higher than 
the homogeneous degree of each term of r. 
Let 

To show that each of the integrals /„ „ is bounded from above, for fixed s such that 
2; ^ we use the estimate of the lower bound of the determinant for (s. A) € Qi {cf. §4.2.2 
and (4.41)), 

\deta.i\>Ci{\Xf + \z\''f. (4.71) 
The principal case 4m — n = 1 is taken care of above see (4.58). For 4m — n > 1 we have 



/ 

Jm. 



dAi 



A" 



deta;!]™ 
= 2Ci 



< 2C- 



dAi 

lAI" 



\XY 



(|A|2 + |0|2)2m 



(lAP + NP) 



2\2m 



dAi + 



lAI" 



where for |Ai| < \z\ we use 
lAp 



< 



N (|AP + |.|Tr» 



< 



dAi 



(|A|2 + |2|2)2m - (|^|2 + |^j^|2)2m-„/2 |^ 

and for |Ai| > \z\ we use the estimate 

|A|" . 1 . 1 



4m— n 



< 



< 



(|A|2 + |2|2)2 - (|Aj^|2 + |Ai|2)2m-n/2 - \Xi\im-n 

Inserting the above estimates into (4.72) yields 

A" 



(4.72) 



(4.73) 



(4.74) 



dAi 



det OL 



;i| 



< 2Cf 



\Am — n 



Cm|^|4m-n-l \^ (4m - n - 1) 



All 
4m — n 



n—Am 



dAi 



(4.75) 



and hence the integral is bounded from above since > and 4m — n — 1 > 0. Thus we 

find that the asymptotic series expansion of r can be integrated, since each term is finite 
for s > Su and args < 7r/2. Furthermore, the A-integral of r-m, which is homogeneous 



B. L. G. Jonsson 



2008-09-04; 11:22 



Wave splitting of Maxwell's equations. 



27 



of order — m, results in b-m + i whieh is homogeneous of order — m + 1. We have hence a 
well defined polyhomogeneous asymptotic expansion of a pseudodifferential operator with 
a parameter of homogeneous degree in {^', s}, the corresponding operator is represented 
in the usual way through an oscillatory integral. 

One can use the residue theorem to evaluate the integrals in terms of the roots of the 
equation detcx-i = 0. This is done for arbitrary Im,n, 4m — n > 1 in Appendix C. 

We have above found an oscillatory integral representation of the desired operator B 
through the A-integral of the symbol expansion of the resolvent. Its principal symbol is 
/ dA r--i, where as usual r--i := cx^^ = (a — A7)~^. One question remains in order 
to associate B with / dA Aj\. It can be reduced to a question of the order of iterated 
integrals. Towards this end we use an alternative representation of the A-integral. We 
note that 

/ dXcx7,^= [ (a;i-J(AR-iAi))-i2(a;i-/AR)(a;i-/(AR+iAi))-i (4.76) 

JxeK ' J \ie[o,oo],XK<T 

where we used the following identity which similar to the (first) resolvent equation: 

(a-i - J(Ar + iAi)) + (a;i - /(Ar - iAi)) 

= (a;i - J(Ar - iAi))"'2(a;i - /AR)(a;i - /(Ar + iAi))"\ (4.77) 

Denote the right-hand side of the above identity w--i{\,S,';s,x). Clearly this identity 
holds also if a-i is replaced with the operator A. By analyticity of the resolvent we can 
choose to integrate along the positive imaginary axis, i.e., Ar = 0. 
Let u be an arbitrary vector in , and consider the two integrals 

Vi{x;s) := I dV ( / dA e'i'--'w,.i{\^'-s,x){Tv){£,',x,)) 

JAr=0,Ai>0 

= ^ei,,[b;o(C', s; x){J^u){e; X3)] (4.78) 

and 

V2{x;s):= [ dX{[ d^ e'«'-"'u;;_i(A, C'; s, x)(.Fu)(e', X3)) 

JXk=0, Xi>0 

= J dXJ^^,lja7,\^',X,s,;x){J^u){^',X3)] (4.79) 

Here we have once again used T to denote the Fourier transform with respect to x' and 
^^i^x:' t*-" denote the inverse Fourier transform from ^' to x' variables. The first integral is 
the standard way of representing the action of principal part of B-o on u. That is B-qu = Vi. 
The second integral V2 is the A-integral of the first term of the parametrix corresponding 
to As,\- We thus have two, possible different, representations of an operator. Below we 
will show that the two representations are equal. For the principal term the problem is 
reduced to showing that the two iterated integrals exist and are equal, e.g., that Vi =¥2- 
We have the following result 

Lemma 2.1. Let u € then for Re{s} > Sr > ^ and args < 7r/2 it follows that 

Vi{-,S) = V2{-,S)^J?. 

This result is shown after the proof of Proposition 2 part 2. We have defined the 
operator B as the oscillatory integral of the A-integral of the symbol representation of the 
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resolvent expansion, and above shown that such an operator exists. The desired sphtting 
matrix is however the A-integral over the oscillatory integral over the resolvent expansion. 
The above lemma shows that the principal term of both these expressions are equal for 
functions on a dense set in the domain. To continue and show that the remaining terms 
in the respective symbol expansions are equal we can once again construct two iterated 
integrals and apply the proof of Lemma 2.1, e.g., the Fubini theorem on this term, and 
since all the assumptions carry over the result remains the same. 

Hence we have shown that the two representations of the splitting matrix indeed are 
equal and can be applied to the wave-splitting procedure below. □ 

4.3 Proof of Proposition 2, part 2 

The symbol of B for the isotropic homogeneous medium case is given in (B.12) (see 
Appendix B) and by counting its powers of it follows that the corresponding operator 
B can be restricted to an unbounded operator on {L^, (•, •)o} with domain and range 
in H". 

To show that this result holds also in general we need to show that the ^'-growth in 
each of the 6-terms is at most linear. To obtain such a result we need good control of 
the shape of the parametrix r oi As,\, which is an asymptotic series of poly- homogeneous 
terms r = r--i + r--2 + r-s + Recall that (see e.g., [32]) r--i := a^^ and 



|„|=i 



r;_2 = -ct7,' [ a;or;_i + ^ [d^^a.,,D2r,_,'\j (4.80) 



r.-m-i = a./ I a-or.-m + a.-ir._rn+i 



|r,|=l 



^ d^^a.,,D2r.,-m+i I , m > 2, (4.81) 

|r,|=2 



where ry G i.e., a multi-index, D^j = \dxj and if ry = (?7i,?72), then DJj. = D!l\D^\ 

Let the linear space of homogeneous polynomials of order n be denoted by hp„. The 
explicit shape of det a-i in Appendix C ensure that det a-i G hp4(s. A, x). Here wc use 
hp„(s. A, a;) to indicate that deta;i is a homogeneous polynomial in s. A, and have 
C°°-coefEcients depending on x. Given the homogeneous polynomials pi G hp„, p2 G hp^ 
we can consider the space of hq^_^ of homogeneous rational fimctions, as elements of 
the form q = pi/p2 and q € hq„_^. We will restrict hq^ even further and require that 
P2 is a power of deta;i. We note two useful properties: Let q G hq„(s. A, x), then 
D^q e hq„(s, A, x) and if qi G hq„, q2 G hq„ then qiq2 G hq„_^„. 

In addition to these two spaces we need two additional spaces. The first is a space 
of block-diagonal matrices hP„, where the two 2x2 blocks have elements which are ho- 
mogeneous polynomials of order n. That is if p G hP„, then p = diag2(Pi, P2), where 
Pi and P2 are 2x2-blocks with each element, {Pm)ij S hp„, for m = 1,2 and i,j = 1,2. 
Clearly for g G hP„(A, ^';a;) and h G hPTO(A, ^';a;) we have hg,gh G \<S'm+n{\^'',x) and 
D^g G hP„(A, ^' , x). The second and final space is W-m and an element h-m is in W-^ 
if it can be written in the form 

4m 

h_„(a,A,e';x) gi':::1f{X,e;x)K-^+''q^_i:^'-\s,X,e;x) (4.82) 

fc=0 neNk 
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where ' G hq^- for a given m, k = 0, ...,4m, and n G Nk and ' G hPj for a 
fixed m, fc = 0, ...,Am, and all n. The matrix if is given in (3.12). Here we have used 
Xm,k,n) index the components of h-m- We require the number of elements 

for each (m, fc)-level to be finite, i.e., \Nk\ < oo where Nk C N. Wc find here the nice 
properties that if h € W-^ then D^h £ W-m and if in addition g € W_„ then we 
have gh G W-n-m- Note that the representation (4.82) is not unique due to that the 
numerator of 9™!^ may contain a power of s. This non- uniqueness of the representation 
will be used constructively in the proof of the lemma below. 
We have the following technical lemma: 

Lemma 2.2. Given the above defined spaces W_m, hp„, and let 

u:= (a, -6), andv.^ (6,6)- (4.83) 

Then, for the homogeneous terms of order —m of the symbol of the parametrix of As,\, 
r--m, we have that r--m € W-m, i-e., 

Am 

r;_„(.,A,r;x) = 9t'^';HA,$';x)A:-+^i':;^'")(.,A,r;x). (4.84) 

fc=0 neNk 

Furthermore, for a fixed m denote the 2x2-block diagonal elements ofg'"^^'"^ by Pin, P^n, 
then Pin = u^win and P2„ = v'^W2n, where each of the elements in the (l,2)-vectors Wjn 
are in hp4^_;^(A, ^';a;) and \Nk\ < oo for each k. 

Proof. A straightforward calculation shows that ct^f^ = r.^i G W_i and r--2 € W_2. 

Furthermore inspection of the leading terms g^^''^'*^'^ and g^^'^''^^ shows that their elements 
can be written as outer products. Indeed, let Pi,P2 denote the 2x2 blocks such that 
=diag2(Pi,P2), then 

Pi = iFu^i and P2 = v'^v^2, where & liP2(Aj C'; x), k ~ 1,2. (4.85) 

Similarly let P3, P4 denote the diagonal 2x2 blocks of 93^'°'°'' then each of these terms are 
of the form u^ws and v-^wa respective, where Wk are (l,2)-vectors with each element in 

The construction of r-^-m-i in (4.81) is a product of finitely many terms, we conse- 
quently find that r.-m-i has at most a finite number of terms. This ensures that the 

\Nk\ < 00 for all k. 

To show that the lemma is valid for an arbitrary m > 3 we make the recursive assump- 
tion that r.-m € W-m and r.-m+i S W-^+i with the desired outer-product structure 
on their respective leading matrices g^^™'^'''^ , 94^-4''''^' fo^' the respective range of j and 
k. 

We now calculate r.^m-i from (4.81) and show that it satisfies the lemma. Towards 
this end wc consider the following three 4x4 matrices c = diag2(ci,C2) G ]iPq{S,']x), 
f = diag2(/i,/2) G hPo{£^';x), and d. The matrix d is defined by d := diag2(M^dii + 
(i[2U,v^ di2 + d22v) where djk are (l,2)-vectors with elements in hpo(^';a;), and u,v are 
the (l,2)-vectors of ^'-elements defined in (4.83). Each term in (4.81) is of the form 

oCi{c+-^dK)hi, where ^1 G W_„ (4.86) 

or of the form 

o.7i-fKh2, where ^2 G W.^+i (4.87) 
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Indeed, the terms containing a-^ and d^a-i with \ri\ = 1 belong to the kind in (4.86), and 
the terms containing a--i, d^a-o, \r]\ = 1 and d^a-i, \r]\ = 2 are of the kind (4.87). 

The lemma follows if we can show that the resulting products of (4.86) and (4.87) 
are elements in W_m-i with the desired outer product-structure on the leading order 
terms. The three terms containing c, d and / respectively are considered separately. For 
the first term we note that chi e W_m- Since a.7i^ e W_i we immediately find that 
a^^chi e W_m_i. The outer-product structure on the leading term survives since the 
leading order term in a~^^ is of the form (4.85). 

Consider the second term containing d. We explicitly write out the leading order 
elements in cx^^^dKhi: 

^ \ " J(-mJm,0,j) , 1 \ " (m,l,n) j^^+i (m,l,n) 

94m ^ 1-4m + ^;;rT 2^ f 4m-l ^ 1-im +' 

j 



1 1 
= + + Po- (4-88) 

Here we let gf™''^' ■* denote the A::th element in hi. 

The first of these terms p„_|_2 = g''^'^'^^ KdK J2k sim °''^^-^'"9^m'''^9-4- The matrices 
in the product all have an outer product structure explicitly given above and from the 
observations that 

KdK = Kdiag2{di,d2)K = diag2(rf2, di) =: d, (4.89) 

and uv'^ — 0, vu^ = Q we find that p.rn+2 = 0- 

To show that s~'^~^p^_^_i +Pq G W_m_i consider 

k,j 

- \^ „(1'0,0) 7 (m,l,«) ^„i+i (m,l,«) (1,0,0) MA.k) ,~im,0,j) (m,0,i) 

— 2_^yi "y4m-l ^ l-im 1-4, + ^ "'94,m ^ l-im 9-4 

n k,j 

= ^Mi„ + ^M2,fc, (4.90) 

n j,k 

where we have once again have used the notation (4.89). Upon multiplying block-diagonal 
matrices with other block diagonal matrices all these with elements which are homoge- 
neous polynomials yield that g^''''''^^ dg'"l^^'''\ 94^'°''^'d34™'ii"^ S hPirn+i- Furthermore, 

(/^m''^9-4 '^^ ^ lic[_4m_4- This suffice for ■s^'^^^Pm+i to be a leading term of an element 
in W-m-i- Similar matrix algebra for a typical term in Pq for a given s-^-order we find 
that each such term fits into a W-m-i element. The remaining issue of the d-containing 
terms is the outer product structure of p^^^. Similarly to the c-terms it is clear that Mi„ 
has the appropriate outer product structure. The term M2jk is a bit more subtle, and we 
need to use the outer-product structure of each of the three matrices. There are two types 
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of terms in g'"^''^'"'^ = P3 +P4, where pg := diag2(pi/,P2-^) with pi,p2 € hp3(A,^';x), and 
P4, := diag2(pnW + u'^pi2,P2iU + v^P22) where pjk are (l,2)-vectors with each element in 
hp2(A,^';x). The pg-term in M2jk yields 



Pil \ fu^dii + dj'^u \ /v^wi 



duv^ wipi 

v'^ di2U^ W2P2 



(4.91) 



and it has the desired outer product structure. To obtain this result we have repeatedly 
used that uv^ = 0, vvF = 0. Similarly for the p^ term we have 



P2lU + v'^P22) V V'^di2+d22v) \ 102 



pfiV + u^p\2 \ ( u^dii + di2'^ \ f v'^wi 

■22 

v?' Pi2U^ dwv^ w\ 

P22V'^ di2U^ W2 



(4.92) 



which once again has the appropriate outer product structure. Wc have hence shown that 
the terms containing d are elements in W-^-i with the desired outer product structure. 
The last kind of terms are these which contain /. These terms are of the form 

' * k 

^ \^ „{rn-lfl,j) T^m-lJm-lfl,]) ^ \^ (m-l,l,n) j^m+2jm-l,l,n) , 

j n 




S' 



+ 9o- (4-93) 



The leading order term is explicitly 

n _ \^„(1.0,0)^^^ (m-l,0j)^m^l(m-l,0j) (1,0,0) . . „ 

j 

Clearly qo := q^^'^+l'^U-f''^^ S hq_4„ and go := 94^'°'°^/94m-4°'^^ ^ liP4„. Observe 
however that qg/ deta-i £ ^'\-im-i that deto:;! = pi}{\,£^' '.x) + s^pi(A, a;) + 
s^P2(A,^';x) where po G hp4, pi S hp2 and p2 G hpg. We note that g^po S hP4m+4 and 
that terms of the form s'^^~"^~^gQPj, j = 1,2 fit nicely as lower order terms in W-m-i- 
Similarly wc can consider the terms of <7q by explicitly calculating the typical -order 
terms and see that each such term fits into W-m-i to finally draw the conclusion that 
ot^^^fKh2 € W-m-i and consequently that r--m-i € W-m-i- The outer product 

structure of the /-terms follows directly from (4.94) and the fact that 94^''^''^^ has the 
appropriate outer product form. 

We can now by a recursion argument draw the conclusion that the lemma is valid for 
all m > 1. □ 

To show that the operator corresponding to (b;-m) maps H"~^ H", it suffices to 
show that (b;-„i)j7c G Sj qi the space of symbols first defined by Hormander. This means 
that we need to show that 

(b;_„),,| < C,,^(l + ICf )i-l''l (4.95) 
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for ri,f3€ N^. From the property of W_„, we know that d!^r-^m € W-m for (s, A) G Qi. 

To show (4.95) for tj = (0,0) we recall that when (s, A) G Qi is deta;i (4.11) elliptic. 
This imply that a typical term of r-m can be bounded as 



Jm 



dA, 



(deta;!)*^ 



<C{\s\' + \en^^, (4.96) 



for all Am — n > 1. For the case 4to — n = 1 the result follows in the sense of a principal 
integral. The symbol b can be written as an expansion of terms which are homogeneous 
in s,^', each constructed by integrating the corresponding r_m-term. The V-m terms 
are of the form (4.84), and for all A; > 1 in (4.84) it follows from each matrix element's 
homogeneity order and (4.96) that the resulting integral is bounded for all values of The 
leading orders A; = have the outer product structure, and wc can therefore apply (4.96) 
to find that the A-integral of these terms grows at most linearly in ^' for large values of 
Applying the derivative dj, on an element in r;_m we find that this reduces the growth 
in ^' with |77|-order for the corresponding b--rn+i symbol, since the integrand is a sum of 
rational functions; it is a sum of polynomials in s,^', A over powers of deta;i both with 
smooth coefficients. The partial derivative is hence a bounded or a more regular function 
in and due to the ellipticity of detcc;! and Re{s} > Sr there exists a A-integrable 
^'-independent function so that we can apply the dominated convergence theorem to show 
that the interchange of A-integral and is allowed. It is clear that the partial derivative 
exists everywhere because that the elements of r-^m are rational functions. We thus find 
that {b--m+i)jk € 'S'l^o and the corresponding operator maps H" to H"~^ as desired. 
This result is valid for any m, and hence we have the same result for B. □ 



4.4 Proof of Lemma 2.1 

To show this result on iterated integrals is a standard application of Fubini's theorem for 
positive integrands, see e.g. [26] . It states that a non-negative measurable function on the 

usual Lcbesgiic-measiirc over the A, ^'-domain has its iterated integrals equal and finite if 
one of them is finite. Wc will apply this to the iterated integrals Vi and V2 with integrand: 

5(A,^';s,a;):=e'«'-^' («;;_! (A, ^',s;a;))i,(.Fw)(0, *,i = l,2 (4.97) 

where 

W-i{X,e,s;x) := (a - 7(Ar - iAi))"'2(a - 7AR)(a - 7(Ar + iAi))"'. (4.98) 

We split the integrand into positive real parts g = gB.+ — gn- + K91+ ~ 9i-)i where ^fe > 
for each k G ii— , /+, /— }. We consider the iterated integral over each of these 

positive functions separately. In order to show that gk for any k G {R-\-,R—,I+,I—} 
is measurable, it is enough to note that w--\ is continuous in both A and ^' and so is 
e'^ , their respective restriction, e.g., the real and non-negative (or any of the other 
combination) is also measurable since they are piecewise continuous. From the assump- 
tions of the lemma we find that u G and it is hence ^'-measurable and consequently 
by trivially extending u to be a function on the product space, it is measurable in both 
A and ^' jointly. We then utilize that products of measurable functions are measurable. 
Consequently we know that the parts of the integrals Vi and V2 corresponding e.g., the 
real positive part exist and are equal. To finish the lemma we note that the operator 
corresponding to b.-m maps to L^. This was shown in the previous section. Hence 

Vi=V2&l?. ' □ 



B. L. G. Jonsson 



2008-09-04; 11:22 



Wave splitting of Maxwell's equations. 



33 



4.5 Proof of Proposition 2, parts 3 and 4 

That B and A commutes in the sense that 

B\iA = AB\^ , (4.99) 

on the set foUows directly form the fact that the resolvent commutes with A due to 
the Fubini-Tonelli theorem and that A commutes in a weak sense with the integral over 
Ai that is used to define B. (See [22] Proposition 2, Part 3.) 

o 2 

To show that B = I one can introduce two projectors defined from the lambda-integral 
over the resolvent, B can be shown to be the difierencc of these two projectors, and the sum 
of the projectors equals the identity. The squaring of the operator is hence the identity 
operator. The key to this proof is to show that the two operators 'P± := ^{I ± B), are 
projectors, it is done by utiUzing the (first) resolvent equation. A proof of the projector 
properties is detailed in [22] Proposition 2, part 3 and 4 for the acoustic case. The 
electromagnetic case follows analogously, and since is somewhat lengthy, we will not repeat 
it here. Once this is known we note that V+ + V- = I, and 'P-'P+ = etc., on an core 

o 2 

set. Consequently, on operator level B = {V+ — V-)^ = V+ + V- = I. □ 

4.6 Proof of Proposition 2, part 5 

o o 

Let 7 ba a scalar, find all (7, £) with non-zero C such that 

Bt^ = ^t^ , (4.100) 
where is a 'vector' of 2 x 2-block matrices of scalar operators. To solve this eigenvalue- 

o 

like problem, we use that B is an involution, that is 

B11S12 + B12B22 = , (4.101) 
+ B12B21 = / . (4.102) 
Writing (4.100) explicitly with 2x2 blocks gives 

Bn£±i+Bi2£±2 = 7.C±i , (4.103) 
B2il^i + hit^i = ^1^2 . (4.104) 
and collecting similar terms yields 

Bi2l^2 = {ll-Bn)l^i, (4.105) 

B2ll^l = hi - B22)t^2 . (4.106) 

Let {Bii +7/) act on (4.105) and use (4.101). Then 



(7'/ - bIi)^^i = {Bu + ll%2t^2 

= Bi2(7/-B22)r^2 

and analogously let B12 act on (4.106) and use (4.102), then 

h2{lI-B22)t^2=h2B2lt^l 

O 2 O I 

= (/-fin)£±i. 



(4.107) 



(4.108) 
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Substituting (4.107) into (4.108) gives after simplification 

(7-727)£=^i =0 . (4.109) 

Thus 7 = ±1 since is assumed to be non-zero. The corresponding eigenvectors are 
obtained by solving the following linear system 



o o 



{±I-Bii)C^i^Bi2C^2 , (4.110) 

{±I-B22)l^2=B2ll^l. (4.111) 

Rewrite (4.102) into 

{±1 - Bn){±I + Bn) ^ B12B21 ■ (4.112) 
Comparison with (4.110) gives that the generalized eigenvectors have the form 

£±=|^±^ + ^"j^± (4.113) 

for arbitrary normalization operators Af^. The condition (B)^ = I imposes, in addition 
to (4.101) and (4.102) the conditions 

B21B11+B22B21 =0, (4.114) 
B21B12 +^22=/. (4.115) 

o 

Hence with use of (4.114) we find that £ is also a solution to (4.111). An alternative 
form of £^ is obtained if we use (4.115) and (4.111) 

which is related to (4.113) by the normahzation {Af"^)' = (±7 - S22)(^=^)" and Af"^ = 
^12 (A/"*)", hence and {C )' differ only by a normalization. 

4.7 Proof of Proposition 2, part 6 

o o 

That B is one-to-one follows directly from the fact that B is an involution. Indeed for any 
7^eC°°(R2,C*) wehave 

eF = 0^F = S0 = 0. (4.117) 

o o 

Hence the null space of B contains only the clement and thus B is one-to-one. That the 

o 

null space is trivial imphes a condition on B21, to see this consider F = (Fi,0) and the 
equation 

BF=|?"^M=0. (4.118) 
\B2iF,J ^ 

00 o 

then, since B is onc-to-onc. 621 is cither trivial or one-to-one, and since the symbol of B21 

o 

is non-trivial, it follows that B21 is non-trivial and hence one-to-one, and thus invertible 
on its range. 
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4.8 Proof of Proposition 2, part 7 

We now extend to a larger domain. This is done by using (S|q)ii and (B|g)2i in the 

o o 

place of Bii and B21 in (4.113), the generalized eigenvector, together with the extension 
of TV"^ to a bounded invertible operator J\f^ : W{m^,C'^) W{R'^,C^) we obtain a 
generalization of to C'^lg : H«(M^,C^) H«^\ where the domain of C'^lg follows 
from the domain of B\^. 



5 Directional decomposition 

We have above collected enough information to proceed and answer the initial question 

about the existence of {£, V}, i.e., docs the decomposition of A exists. Most of the proof 
are done for the set C°°(M^, C^), but in the end we extend the results to the general case. 

With the definition of the splitting matrix in Proposition 2, in particular the commuta- 
tion between the splitting matrix and the electromagnetic system's matrix (see Proposition 
2, part 3), wo obtain the decomposition by the following proposition. 

Proposition 3. The equation 

AjC\^ = £|iV , 

has a solution where the columns of are the generalized eigenvectors ^^\q to for 
q = 1,3 and where V is a block diagonal matrix with the elements : H^(M^,C^) — > 
H^(M'^,C^), representing a generalization of the vertical wave number and 

= {J^\,)-\{B\,hi)-' {A2i{±I - (e|3)n) + ^22(^13)21) AAI3 , (5.1) 
where J\f\g = J\f+\g = J\f-\g, 9 = 1,3. 

Remark 3.1. If one considers the equation AL'^ = £^<S^ with the particular normal- 

o 0-1 

ization M = B21 and let 

^ i±I + Bn)B2i , (5.2) 

o . 

then, upon eliminating S^, one finds that establishing the decomposition is equivalent to 
solving the equation 

OiO 0_|_ OiO o o 

Z^A2iZ^ +Z^A22- AiiZ^ - Ai2 = , (5.3) 
i.e., an algebraic Riccati operator equation. As C solves the decomposition problem 

° -ir- 

we have the fact that Z^ solves the associated algebraic Riccati operator equation. The 
map Z"^ is denoted the impedance mapping. Note that one can obtain a corresponding 
admittance mapping, {Z^)~^ , that solves the algebraic Riccati operator equation that is 
obtained by operating with {Z'^)~^ on both sides of (5.3). 



5.1 Proof of Proposition 3 

To show that C.'^\q decomposes^, webegin with the proof on C°°(M^,C^). By Proposition 
2, part 3 and 5 we have 

000. 000. 00. 

BAC^ = ABC^ = ±AC^ . (5.4) 
Let H± := AC^, then by (5.4) 

BH^ = ±n^ , (5.5) 
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but from Proposition 2, part 5 we know that for some arbitrary normalization operator, 
A^^, we have 

H^=(^^ + ^''^M^ = l^S^, (5.6) 

I o o o 

by the choice of normaUzation operator M = Af S"^, for some particular <S*. From the 
definition of H^, we get 

aI^ =n^ = l^S^ . (5.7) 

Hence C are generalized eigenvectors of A with generalized eigenvectors S of 2 x 2 

o 

matrices of scalar operators. To obtain an explicit expression for iS*, consider (5.7) 
explicitly 

(in(±/ + ^ii) + Ai2B2i)J<f^ = (±/ + Bn)J^^S^ , (5.8) 

{A2i{±I + Bii) + A22B2i)M^ = B2iM^S^ . (5.9) 

From (5.9) wc obtain that the range of the left and right hand sides have to agree, thus 

o 

the left hand side is within the range of 621 find hence the inverse is defined on this range 
and thus 

5± = (^±)-i(S2i)-i(i2i(±/ + Bn) + A22B2iW^ , (5.10) 
is well defined and it also is the generalized eigenvalue to A. To see that (5.8) gives the 

o o 

same result, we use two of the equations implied by the commutation of A and B, viz. 

AuBu + A12B21 = BuAn + B12A21 , (5.11) 
A2iBn + A22B21 = Bsiiii + B22A21 . (5.12) 

We rewrite (5.8) and apply (5.11) to obtain 

(±iii + AuBii + Ai2B2iW^ = ((±/ + ^ii)iii + Bi2A2i)M^ 

O O I O I 

= (±/ + 6ii)AA±5± . 
Applying ±7 — Bn to both sides and using (4.101) and (4.102) gives 



(5.13) 



((/ - B^,)An + {±1 - Bn)Bi2A2iW^ = {I - B,,)Af^S^ ^ (5.14) 
Bi2{±A2i + B2iAn + B22A2i)Af^ = BuB2iJ^^S^ , (5.15) 
and using (5.12) gives 

Bi2(i2i(±/ + Bn) + A22B2i)Af^ = B^kiAf^S^ , (5.16) 

and hence an equation for S that is equivalent to (5.9). Thus we have shown that the 
two expressions (5.8) and (5.9) are equivalent and that (5.10) is the solution to both. 
Before we extend (5.7) to the general one, we introduce the matrix operators 

V=P" P] andr=f^t^- -^o+M^, (5.17) 
\0S- \ B21 S21 / ^ ' 
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where we have made the choice Af = Af~^ = N . With the introduced notation, (5.7) 
becomes 

At = tv , (5.18) 
and furthermore, we may rewrite this equation into 

AL=L\^V, (5.19) 

where by the definition of we know that is well defined, the choice of C\-^ follows 
since : To extend the domain to a larger set, let {Fn}^^-^ C C°°{R'^,C^) 

be a Cauchy sequence. Then for fixed r] > and large enough n we have 

\\ACFn-ACLF4o<V, (5.20) 



as long as lim„^oo Fn = Fq G D(£|g) = H^. Let 



3 ' 

AC\^Fo = G. (5.21) 



Subtracting (5.21) from (5.19) and using (5.20) we obtain that for large enough n 

V > \\ACF„ - A i^ollo = II VFr, - G\\o , (5.22) 

o 

thus the hmit of the right hand side exists, that is there exists an extension of V. 

o 

Due to the requirement of equal domains of the extensions of V and ^ ^Cjg we find 
that 

G=C\,V\sFo, (5.23) 
where the elements of Vjg, have the form 

5^ = (^|i)-^((S|,)2l)-^(^2l(±/ + (S|3)n) + ^22(S|3)2l) Af\, (5.24) 

and 5± : H3(M^C2) ^ H1(IR2,C2). Thus on we have obtained 

^'C|3 = £|iV|3 . (5.25) 



6 Discussion of the result 

By applying functional analysis to the problem of decomposition of the wave field for 

the electromagnetic system's matrix we have extended the wave-splitting procedure to an 
anisotropic media whose properties vary with all three spatial coordinates. The result 
extends beyond the up/down symmetric case. The analysis of the spectrum shows that 
a strip around the imaginary axis is in the resolvent set. We define a resolvent integral 
whose contour lies in this strip. Due to the explicit form of the systems matrix, we do 
not have a full spectral resolution of the operator. StiU, the resolvent integral over a 
path in the resolvent strip is shown to be well defined by applying the elliptic theory of 
pseudodifferential operators with parameters. 

Using this resolvent integral we define a splitting matrix. This matrix has the fea- 
ture that one can construct its generalized eigenvectors of operators corresponding to the 
(generalized) eigenvalues ±1. We have above shown that the splitting matrix commutes 
with the electromagnetic system's matrix. One consequence of this 'commutation' of the 
operators is that the generalized eigenvectors of the splitting matrix also are general- 
ized eigenvectors to the electromagnetic system's matrix. The corresponding generalized 
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eigenvalue to the system's matrix (a 2 x 2 matrix operator) is the key ingredient in the 
definition of one-way equations for the electromagnetic case. This 'eigenvalue' is the elec- 
tromagnetic generalization of the vertical wave number obtained in the linear acoustic 
case. One of the features of this procedure of decomposition is that we have constructed 
the composition operator without having to invert any of the elements of the splitting 
matrix: the construction relies on the fact that the splitting matrix is an involution than 
in the corresponding acoustic case. However, this result can also be carried over to the 
acoustic case. The removal of the inverse of an element in the splitting matrix from the 
splitting process, is not complete. It remains in the one-way equation obtained after the 
splitting, even though we have been able to remove it from the composition operator. 
The generalized eigenvectors to the splitting matrix is used to generate the composition 
matrix that decomposes the electromagnetic system's matrix. 

The traditional approach to the decomposition of the system's matrix gives an alge- 
braic Riccati operator equation. The splitting matrix construction of the decomposition 
gives us a family of solutions to this operator equation in terms of the elements of an 
integral over the resolvent of the electromagnetic system's matrix. 

Once we have obtained the wave decomposition, we can proceed and use the one-way 
representation to study direct problems or by applying the generalized Bremmer coupling 
series to study direct and inverse scattering problems. 



A Derivation of A 

Given the normalized Maxwell equations (2.5) in a medium with the constitutive relations 
(2.3), we have Eq. (2.6) with vertical components (2.7) and transverse coordinates (2.9). 
We replace the explicit appearance of sH^ and sE^ with (2.7) and obtain 



(A.1) 



SVc^pHp - A<a3M33^(V X E)^ + (V X E)^ = K°- fXaSfJ-siK^ , 

-seocf}Ef3 - eas^ssiy X H)3 + {V x H)^ = Jl- f^asf^ss^i ■ 
where 

SafS = ^afj — easesa^es/J , I'afj = fJ-aff — IJ-aSfJ-si IJ-ais ■ (-^-2) 
The transverse components of V x E and V x H are explicitly 



(V X E)i = d2E3 - d^Ei , (V X H)i = d2Hs - 53/^2 , 

(V X E)2 = dsEi - diEs , (V X H)2 = d^H^ - diH^ . 



(A.3) 



We replace the vertical components (£^3,^3) in (A.3) with the identity (2.7) to find 

(V X E)i = -d3E2 - 92e3~3 esa^a + 826^^ {diH2 ~ ^ai^i) " s-^d2e^^^JI , 
(V X E)2 = dsE, + die-^hs^E^ - s'^ d^e-^ {d^H2 - Sai/i) + s-^d^e^^J^ , 

(V X iJ)i = -d3H2 - 52^33^3/3-^/3 - 821133 {d^E2 - d2Ei) + 821133 Kl , 

(V X H)2 = d3Hi + ^x^i3^^x30H0 + s-^^x^i3^{^xE2 - d2E^) - s-'^difisixi . 

These transverse components of the curl is substituted back into (A.l). Collecting similar 
terms gives 

SV213H13 - ll23lJ'33{dlE2 - d2Ei) + dsEi + die^3€3aEa 

- s-'die33\d,H2 - d2H,) = iq - fi23fi3iKI - s-^^ieg/jg^ , (A.4) 
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+ s-^d2e^s\diH2 - d2Hi) = - MiaMsa'^l + s-'aaegg' J| , (A.5) 

sexpEp + eise^^{d\H2 - (?2H"i) + 93-^2 + ^aMsaVs/S-ff/s 

+ s-^^2^Xsii^lE2 - d2Ei) = -Jt + eiaeia^ J| + s-^SaMia'-^l - (A.6) 



- SS20E0 - e23e^^{diH2 - d2Hi) + d^Hi + dni^^^ ii^jsH p 

+ s-'SiMas' {diE2 - d2Ei) = J| - easeis'^l + s-^d^ix^^Kl . (A.7) 

Separating the derivatives in the vertical direction from the remaining terms gives 

{Ids + A)F = N , (A.8) 

in which the elements of the electromagnetic field matrix, F, arc given by 

Fi := El , F2 := -E2 and F3 := H2 , F4 := Hi . (A.9) 

We write the electromagnetic system's matrix, A, as a matrix of 2x2 block matrices 



A:= 

where each block-matrix is given by 

All 



All A12 
A21 A22 



(A.IO) 



^13^2 A*l3i9l 



A12 

A21 
A22 



:= s 



V22 V21 
£11 — £12 

—£21 £22 

92^32 92^31 
d\[i.Z2 difl3i 



— S 



dif^i — 9ie32 
—92631 92632 

die^.^di -die^^d2 
d2c^idi d2e^zd2 

-1 ( 92^33^92 92^33^ 9i 

9iMi3^92 din^^di 

eisdi —£1382 

-€23dl £23 92 



■^33 ' 



(A.11) 



A'ss + ^33 



and the elements of the source terms 

K2 - ^^23^^33Kl - s^^die^^JI , 
Kl - Mi3M33^^s^3 + 5^^92633^3^ , 
+ 613633^3° + 5-193^33^^3 ' 
J2 - £23633^3" + . 



Ni 
N2 

N4 



(A.12) 



Note that the adjoint of the 2x2-block matrices with respect to the L^-inner product 
satisfy the following relations 



for self- adjoint e, /x. 



-^22 , (^12 (s))* = ^12(5) , (^21 (s))* = ^21 (s) , 



(A.13) 
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B The isotropic homogeneous case 

We consider the normalized Maxwell equations (2.5) for isotropic homogeneous media, 
i.e., we assume the constitutive relations 

B{x,s) ^ fiH{x,s) , 

D{x,s) = eE{x,s) , ^ ' ' 

where fi, the permeability, and e, the permittivity, arc both real valued scalars and inde- 
pendent of space and time-Laplace parameter s. The electromagnetic wave field satisfy 
Maxwell equations 

^ B.2 
-esE + VxH=r . ^ ^ 

The electromagnetic system's matrix, derived in Appendix A, cf. (A. 11) reduce in the 
homogeneous isotropic case to 



11 = , 

T „-l,-l ( 

8281 8282 

8282 -c 

8182 8181 



A X -1 -1 / 8181 8182 

Ai2 = sfil - s ^ ' 



A21 = sel — s jjL 
^22 = , 



, (B.3) 

1 _i / 52^2 -8281 



where I is the 2x2 identity matrix. Applying the fourier transform in transverse space, 
and using that the coefHcients are constant, gives the spectrum as the set of A(^', s) such 
that det(a(^', s) - A) = or 

(s^e/x+ICf -A2)2 = ^A = ±V^VTW- (B.4) 

Note that Re{s} > imply |Ar| > 0. Hence for Re{s} > the spectrum separates into 
two parts. The inverse of a;i(C'j s. A) := a(^', s) — A/, which we denote with r;_i(^', s. A) 
is 

r,_i(r, .s, A) = s, -\){s^e,ji + - X^Y' . (B.5) 

The splitting matrix defined as, cf. (4.2) 

6;o = ^ / a7^^ , (B.6) 

has two parts, one is proportional to 

1 / dA ^ 1 .0 7^ 

TTi 7,^=0 + If |2 - A2) ^/iVTTP ■ 

and the second part is 

as a principal value. With the observation that (ai2(f , s))"-^ = a2i{£,' , s)/{s'^eiJ + |f P) 
and (B.7)-(B.8) we obtain 

b;o(^',s)=f °i i) , (B.9) 
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where 

Note that b-o has the expected property: = / cf. Proposition 2, part 4. 
As the material is homogeneous we have 

B = jr-i6.o. (B.ll) 

The symbol corresponding to the generalized eigenvector, Proposition 2, part 5, is 

= (^^/i) n± , (B.12) 

for some normalization n^. The generalized eigenvalue of A corresponding to £^ has the 
symbol representation cf. Proposition 3 

= ±n~^za2in = ±n~^ ^==^===012021/1 



= ±x/sV+le?/ • (B.13) 

Thus the generalized eigenvahie problem reduces to an eigenvalue problem for the isotropic 
case, i.e., reduces to diagonal matrices. If we consider the wave-splitting problem by 
earlier developed techniques see e.g. [8, 30] we find 

= (B.U) 

Hence the general procedure described in this paper agrees with the earlier wave-splitting 
methods available for the homogeneous (and layered homogeneous cases. 

C Two tools for the proof of Proposition 2 

C.l The determinant of a.i 

The determinant of the symbol of A — XI, is 

det «;! = (A^ - 2ie33^Re {esa}^aX - ea/3^a^/3) (A^ - 2iij,^^Re {/j.^^} ^^A 

- + det sdetu + [A^(£i2Z^i2 - £222^11 - £llt'22 + £2lt'2l) 

+£a/3CaC/3M33^ det + ValB^a^^esi det € 
-2(e33Ai33)"^R.e{(e3. X ea.)3(M3- x iy,3-)3^a£.l3} 
-2iA (Re {{{11.3 X K:)3 X e-.ah^a} + Re {((e.3 X e.:)3 X !/:a)3Ca})] , (C.l) 



where 
and 



(e3. X £a.)3 = e3i£a2 — e32£al (C-2) 



((/i.3 X U..,)3 X £:a)3 = /Ui3l'2l£2a " M232^1l£2a " Ml32^22£la + M232^12£la • (C.3) 

The similar terms are defined analogously. Notice that (C.l) for Ar = this is a second 
order polynomial in with real coefficients. 



B. L. G. Jonsson 



2008-09-04; 11:22 



Wave splitting of Maxwell's equations. 



42 



C.2 Partial result for the symbol of the splitting matrix in the 
anisotropic case 

To obtain the symbol of the sphtting matrix, integration of the type Im,n cf- Section 4.2.4 
Eq. (4.69), is to be evaluated. In the polyhomogeneous expansion it is clear that the case 
4m = n + 1 exists only for m = l.n = 3 duo to the homogeneous decreasing degree of 
the polyhomogeneous expansion of r. This is the principal integral and it is calculated 
in Proposition 2, part 1. For the remaining terms, 4rn > n + 1, we use homogeneity of 
det a-i to obtain. 



lim 

A I — * oo 



(det cy.-i 

that allows us to use the residue theorem. Thus 



< lim 2|A|-^ 

|A|— i-oo 



(C.4) 



27r J_ 



A" 



<iAi / , , ^ 
^ (deta;i) 



= —Res 



A" 



(det a- 



Xt \ - Res 



(deta;i)™' ^ 



where the roots of the fourth order polynomial det a-i are denoted by Aj^ and Aj . In the 
evaluation of the integral we have to consider the case when \ f = A J . 

To find the residue at Xf we first consider the case that Af ^ A2 and choose A such 



that 



Use the identity 



A- 


-A+ 


< 1 and 


A- 


-A+ 




-A^ 


A^ 


-Ar 



< 1 . 



(C.5) 



(m- 1)! 
(1-2/)'" 



dt/'^ 



2^ — 7^ — y 



(C.6) 



3=0 



valid for \y\ < 1 and m G {1, 2, 3, . . .}, together with the binomial theorem to rewrite the 
integrand of Im,n into the following form 



A" 



(det a;!)™ 



p=0 



{X- xt)p{xt)''-p 



= E 



(A - A+)™(A - A+)'"(A - Ar)™(A - A2-)'" 
n\ (A+)"-P 1 ^ (_i)ji+i2+j3 



^„ vw (A - xtr iim - ly-r ,^,,^-=0 iiii^bai 

{m. + ji — 1)!(to + j2 — l)!(m + ,73 — 1)! 



(C.7) 



—{x-x+r+^' 



+j2+j3 



(A+ - A+)^i+™(A+ - Ar)^=+™(A+ - A2-)J3H 
The residue at A^" is the coefficient of the sum such that ji + j2 + ja = m — p— 1. Thus 



Res 



A" 



(dettt;!)*^ 



p=0 



(A^)' 



pj ((m- 1)!)3 



E 



ji +32 +33 =m-p- 1 



(^^iyi+32+j3 

ii!i2!i3! 



{m+ji - l)!(m + j2 - 1)!(to + j3 - 1)! 
(A+ - A+)Ji+'"(A+ - X^y^+'^iXt - A^)J3+m • 

An analogous result is obtained for the root A2 • Observe that the term (A^'" — A^)^-'^"™, 
is not bounded, but from Proposition 2, part 1 we know that the integral is bounded, and 
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hence this can be removed by ehminating common factors in the sum of the two residues, 
similarly to (4.62) and (4.63). 

For the case of = we obtain that the two residue collapse to one and becomes 

'''n(deta;i)'"'^M ^oW((m-l)!)2 

^ V- {-iy^+^Hm+,n-l)\{m + j2-iy. 

Jib2!(A+ - Ar)-+"(A+ - A,-)-+™ • 

ji>0 

Thus given the roots of the polynomial detctji = 0, we obtain the integral for each 
m, n. Upon substituting the integral in the asymptotic series for the parametrix we obtain 
the symbol. 
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